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Introduction.^ 

The present memoir was originally intended to deal only with the theory of diffusion, 
which still remains its chief subject. During the course of the work, however, it 
became clear that the theory of viscosity and thermal conduction could also be 
incorporated by a slight extension of the analysis. This has been done, and the 
paper now affords an account of all these three ''' ordinary '' mean -free-path phenomena 
of a composite gas. 

The treatment of viscosity and conduction is brief, partly because the theory 
for a composite gas is so much more complex and less important than that for 

* See Note F, p, 197, 
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116 DR. S. CHAPMAN ON THE KINETIC THEORY OF A COMPOSITE 

a simple gas.^ For the same reason the results are not carried to a higher degree 
of approximation than that attained, in regard to the same phenomena, in an earlier 
memoir, t The present method, however, enables the approximation to be carried 
to any degree of accuracy, which was not formerly possible. Also a certain mistake 
in the previous investigation of the conductivity of a composite gas is indicated and 
corrected {cf. § 18). 

By its very nature, the problem of diffusion requires the consideration of molecules 
of two kinds, a complication which is unnecessary in the construction of a theory of 
viscosity and conduction. Perhaps this circumstance largely explains the greater pro- 
gress which had hitherto been made in the latter theory, as compared with the theory 
of diffusion. Until recently the only accurate expression which had been obtained 
for the coefficient of diffusion D^^ was that deduced by Maxwell in his second great 
memoir J on the dynamical theory ; it referred exclusively to a gas whose molecules 
inter-act according to the inverse fifth-power law of force. The best available value 
of Di2 for molecules of other kinds was due to Langevin,§ but the formula, unlike 
Maxwell's, was only approximate. It was determined on the assumption that the 
distribution of velocities in each group of molecules, relative to the mean velocity 
of the group, was according to Maxwell's well-known law appropriate to a gas 
in the uniform steady state. The amount of error (if any) introduced by this 
assumption was unknown. In the present paper the true law of distribution is 
determined, and an exact expression is obtained for D12 which is applicable to the 
most general case of a composite monatomic gas. It is found on comparison that 
the error of the above approximate formula is as great as 13 per cent, in extreme 
cases (§13 (e)). 

A particular case of Langevin's formula, relating to rigid elastic spherical 
molecules, had previously been deduced by Stefan in 1871. The theories of 
MaxwellJI STEFAN,ir BoLTZMANN,^'^' and Langevin, and my own earlier theory, 
all agreed in predicting no change in D^g with the relative proportion of the two sets 
of molecules. Another theory, originated by MEYEuft, asserted that there would 
be a large variation in J) 12 as the proportion of either component varied from to 1. 

^ This has been dealt with in detail in my recent memoir, * Phil. Trans.,' A, vol. 216, pp. 279-348, 1915. 

t 'Phil. Trans./ A, vol. 211, pp. 433-483, 1911. 

X Maxwell, ' Collected Works,' ii., p. 27. His formula for D12 is a special case of the general 
result (13'03) of this paper. 

§ Langevin, ' Ann. de Chimie et de Physique,' (8), v., 245 (1905); cf, also Enskog, 'Phys. Zeit./ xii., 
533 (1911). The same result was independently discovered by the present writer, * Phil. Trans.,' A 
vol. 211, p. 499 (1911). 

II Maxwell, ' Collected Works,' i., p. 392 ; ii., p. 57, p. 345. 

U Stefan, ' Wien. Sitzb.,' 63, (2), p. 63, 1871 ; 65, p. 323, 1872. 
*^ Boltzmann, ' Wien. Sitzb.,' 66, p. 324, 1872; 78, p. 733, 1878; 86, p. 63, 1882; 88, p. 835, 1883 ; 
also ' Vorlesungen,' i., p. 96. 
tt Meyer, * Kinetic Theory of Gases,' p. 255 (English ed.) ; also Gross, * Wied. Ann.,' 40, p. 424, 1890. 
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The experimental data which were determined in order to test the point confirmed 
neither conclusion ; a variation in Dj2 was observed, but it was so much smaller than 
that suggested by Meyer's formula that the experimental values of D12 were in closer 
numerical agreement with a constant expression than with Meyer's highly variable 
result. Eecently Kuenen^ has modified Meyer's theory by taking into account 
the tendency of a molecule to persist in its original direction after collision, a 
consideration the importance of which was brought into prominence by JEANsf 
in connection with the elastic-sphere theory of viscosity. As a result of this revision 
Meyer's theory is brought much more closely into accord with experiment (§13 (i) ), 
though the discrepancies still exceed those between the observed values of D12 and 
a suitable constant mean. 

In § 13 (i) the observations of D12 relating to the above variations are compared 
also with the results of the present theory, which affords a formula similar to that 
of Langevin but multiplied by a correction factor which is a function of the relative 
proportions of the component gases. The agreement with experiment, while not 
exact, is perhaps as close as the degree of experimental errors, and the uncertainty 
as to the best molecular model, entitle us to expect. 

The present theory of diffusion is compared also, in two particular cases, with 
exact results obtained theoretically by other writers. These cases are (a) that of 
a gas formed of two sets of molecules which are identical in their dynamical 
properties (as in the problem of the self diffusion of a gas), and (&) that of a gas 
in which the molecules of one kind are infinitesimal in size and mass compared with 
those of the other kind (as in electronic diffusion in gases). These problems have 
been dealt with respectively by PidduokJ and Lorentz§ ; the corresponding special 
forms of the general expressions here . deduced are in satisfactory accordance with 
their results (§ 13 (e), (/)). 

Certain other phenomena connected with diffusion are also discussed which do not 
seem to have been investigated in any detail hitherto. In §§ 10, 14, 15 it is shown 
that diffusion may be produced by a gradient of pressure or temperature independently 
of the presence of a concentration gradient or of external forces; in §§14, 15 the 
amount of this effect is considered, with numerical illustrations relating to particular 
pairs of gases. Conversely, the absence of diffusion in a composite gas in which the 
temperature is non-uniform (the pressure being uniform and there being no external 
forces) implies the existence of a corresponding variation in the relative concentration 
of the two gases. This latter variation (c/. § 16) appears to be so large in amount 

* KuENEN, ' Supp. No. 8 to the Communications from the Ley den Physical Laboratory/ January, 1913 ; 
also * Amsterdam Acad. Proc.,' 16, p. 1162, 1914. 

t Jeans, ' Dynamical Theory of Gases ' (2nd ed.), pp. ^76, 292, 328. 

t PmnucK, *Proc. Lond. Math. Soc' (2), 15, p. 89, 1915. 

§ LoRENTZ, * Archives N^erlandaises,' 10, p. 336, 1905; 'Theory of Electrons,' p. 268. More general 
results obtained by Jeans, using Lorentz's method, are given in his 'Dynamical Theory of Gases,' 
2nd ed., §§ 314, et seq. Of. also PmnuCK, he, cit, p. 112. 
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as to suggest a doubt whether some compensating influence has not been overlooked 
in the theory, and it is desirable that the matter should be put to the test of 
experiment (see Note A, p. 196). In §§12, 19, it is shown that diffusion is neces- 
sarily accompanied by a transfer of thermal energy, and a new physical constant, 
the *^ specific energy of diffusion,'' is introduced and discussed. 

The method used to determine the velocity-distribution function is similar to that 
published in my recent memoir on a simple gas ; the details of the work are, of 
course, more complicated in the present case. The formulation of the equations 
of diffusion and energy for a composite gas, executed in §§10 and 12, embodies 
certain features which seem to be novel. 

We may here remark also upon some by-products of the analysis which suggest 
interesting developments in the field of pure mathematics. The comparison of 
LoRENTz's solution of the problem of electronic diffusion with my own has led to 
expressions for ir and sin ttx of an altogether new form. Lorentz used Boltzmann's 
integral equation for the velocity- distribution function, and obtained a solution in 
finite terms involving tt ; the solution arrived at in this paper is determined by the 
use of the aggregate of the equations of transfer (§ 2), which is really equivalent 
to Boltzmann's equation. The result is expressed, however, in terms of the quotient 
of a symmetrical infinite determinant by its principal minor, and formulse of this 
kind are hence found for tt (and also for sin irx). The elements of the determinant 
are expressible simply, in terms of gamma functions. A further study of the subject 
from the analytical point of view would probably be fruitful in results of interest and 
importance (see Note B, p. 196). 

I hope later to apply the present methods to the examination of the problems 
offered by rarefied gases. So far, however, as concerns the mean-free-path phenomena 
in monatomic gases under normal conditions, the investigation imperfectly attempted 
in my memoir of 1911 is completed by this and the second paper already referred to 
.('Phil. Trans.,' A, vol. 216). 

It is a pleasure here to make grateful acknowledgment of my indebtedness 
to Sir Joseph Larmor throughout these investigations, which were started under 
his influence, and would hardly have been carried to this stage but for the inspiration 
afforded by his continued encouragement and interest. 

§ 1. Analysis of the Dynamical State of a Composite Gas. 

(a) Notation,'^ 

The gas considered in this paper is one composed of molecules of two kinds, each 

having the property of spherical symmetry (or, in brief, each being "monatomic"). 

■^ In numbering the equations I have adopted the decimal method introduced by Peano. The 
number to the left of the decimal point is the number of the section, and within any section the numbers 
to the right, if read as decimals, are in numerical order. With this method it is possible, by the 
introduction of a third or even fourth figure, to number equations inserted between others already 
numbered, without having to alter the references to all succeeding equations. 
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The density of the gas is supposed to be such that the mean free path of a molecule 
is large compared with the distance at which molecules appreciably affect one 
another's trajectories : this we express concisely by terming the gas *' nearly 
perfect." 

Similar quantities relating to the two groups of molecules will be represented 
by similar symbols, with distinguishing suffixes 1, 2 ; it is conyenient to adopt the 
convention that the first gas is that which has the gi^eater molecular mass. The 
molecular masses will be denoted by mi, mg, while the notation for various other 
characteristics of the gas at (x^ t/, z, t)^ is explained by the following list : 

vi, V2 ^ the number of molecules of the first and second kinds 
per unit volume. 

X3, X2 = the proportion of molecules of each kind at {x^ y, z^ t), 

Pi, p2 = the densities of the constituent gases. 

Pi, P2 = the external forces (in vector notation) acting on each 

molecule mi, mg. 

(Xi, Yi, Zi), (X2, Yg, 7i^ = the same forces in Cartesian notation. 

Ci, C2 or {ui, Vu ^(h)y (^2? '^2> '^s) = the mean velocities of the two groups of molecules in 

vector or Cartesian notation. 

We define further quantities of the same nature, in terms of the above, as 
follows : — 

(rOl) vq = V1 + V2, so that Xi = vi/{vi-\-v2) = ^iAo> ^2 = ^2/(^1 + ^2) = ^2/^0? X1 + X2 = 1. 
(1'02) 2XV= X1--X2, X12 = Xi/Xg, X21 ~ X2/X1, so that Xi2X2i = 1. 

(1'03) rriQ = Ximi + Xgrna, m'o = Xi (mi-mo) = -Xg (ma-mo) === XiX2(mi— ma). 
Now we have 

(1'04) >i = Pimi, p2 = 1^2^25 

so that, by (I'Ol) and (l'03), 

(1'05) Po = Pi + P2 = ^1^1 + ^2^2 = ^0^0, /o = T~ = ^o(^l~^2)- 

^1^2 

Also, in vector notation P, c, we shall write 

( 1 '06 ) Pa ^ XiPi + X2P2? Pa = Xi ( Pi Pq ) = — • X2 { X 2~~ — - 1 n ) ~ X1X2 ( Pi P2 )• 

, \ mo / \ mo / \^?^o ^^0 

(1-07) Co ~ X1C1 + X2C2, c'o = Xi(ci~Co) = -X2(c2"-Co) = X1X2 (ci-Cg). 



^ 



I.e., at the point (x, y, z) and at time t. 



120 DE. S. CHAPMAN ON THE KINETIC THEORY OF A COMPOSITE 

The corresponding component equations, involving X, Y, Z or u^ v, w with 
appropriate suffixes, are similar and will not be written down here. 

By inverting the above equations we obtain the following expressions for the 
original in terms of the derived quantities : — 

(1-08) ^^=.^^+!^o^ p^ = ^ip^+^^ c, = Oo+r^ 

Ai wIq Xi Xi 

(r09) rri.^m,^'^, p^^^^p^^^lo^ c.^c,^^^ 

Xo mn Xo 



2 



(b) The Interpretation of the Derived Quantities. 

Corresponding to (I'OS) and (1*09), the motion of the gas can be analysed into (a) 
a steady motion of the composite gas as a whole with velocity Cq, together with (6) a 
motion of interdiffusion in which the mean velocities of the two streams are 
respectively c^o/Xi and — cyXg. In this latter motion equal numbers of molecules are 
transferred per unit time in each direction, the number (per unit area normal to the 
direction of the vector c^o) being vqC^q^ since by (roi) 

(1*10) n^yXi = V2C^o/^2 = ^0<^V 

The momentum of the common motion (a) is clearly pqCq per unit volume, while that 
of the motion of interdiffusion is i^o(^i"~^2)^^o ^r, by (l'05), pqC^qI in general this is 
not zero, owing to the inequality of mass of the molecules, although the diffusing 
streams convey equal numbers of them in any interval of time. This analysis of 
momentum corresponds to the equation 

(111) piCi + P2C2 = i/imjCi + v2m^C2 = po^o + P W 

The equations (1*06) differ from the others by involving the molecular masses as 
well as vi and vq. This resolution of the forces Pj and P^ may be considered as 

follows : the first terms (c/. 1*08 and 1*09), viz., — Po on mj and — Po on m^, represent 

mo mo 

forces which will impart a common acceleration Po/mo to each group of molecules (we 
may regard this variation as affecting their common velocity of streaming, Cq) ; the 
remaining components, PVXi on wii and — PVXg on mg, when summed up over the 
Pi, 1^2 molecules of the corresponding groups, afford equal and opposite total forces 
i/qP^o, vqP%. In connection with this we may remember that two such interdiffusing 
groups of molecules as we have considered will exert equal and opposite actions on 
one another, and that equal and opposite forces must be applied to the two groups if 
their motion of interdiffusion is to be maintained, or modified without imparting any 
common velocity to them. 
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As regards X'o, this Is useful for the purpose of imparting a symmetrical form to 
some of our equations ; if D denotes a differential operator of any kind with respect 
to (x, y^ Zyt)^ since X1+X2 = 1 we have, by (l'02), 

(112) DXi = -DXs = ID (X1-X2) = DxV 

(c) The ^' Peculiar'' Motions of the Molecules. 

So far we have been concerned with the mean velocities of the constituents of the 
gas, without considering the actual motions of individual molecules. We shall denote 
the velocities of typical molecules of the two kinds by {0)1,(0)2 or {{u)i, {v)i, {w)y} , 
{(te)2, {v)2y (^)2}> when referred to the co-ordinate axes, or, when referred to axes 
moving with the velocity Cq appropriate to the point and time in question, by Cj, C2 
or (Ui, vi. Wi), (U2, V2, W2). Thus 

(113) Ci = (c)i~Co, €2 = {0)2-00. 

The velocities Ci, Cg will be called the " peculiar " velocities of the molecules. 

We have no means of determining the individual values of Ci, Cg for the molecules 
near {x, y^ %, i), but for a given state of the gas, as specified by its composition, mean 
motions, pressures, and temperature (these are expressible in terms of mean values of 
functions of C), there will be a certain frequency law, or function representing the 
distribution of various values of the velocity among the molecules. The determination 
of this velocity -distribution function is fundamental in the method of this paper. It 
will clearly involve C or its components as independent variables, together with 
certain parameters {e.g., pressures or mean velocities) which are dependent on 

The mean value of any function of the molecular velocities will be denoted by 
placing a bar over the expression representing the function. Thus, for instance 

{of. ro7), 



(114) (c)i = Ci, {0)2 = 02, Ci = (c)i-Co = Ci-Co = c'o/Xi, C2 = {0)2-00 = 02-Co=-c^o/K 

It is convenient at this stage to modify the meaning of our symbols Ci, Cg, which 
have so far represented vector quantities ; henceforward they will denote not the 
vectors themselves, but their amplitudes. These, of course, are essentially positive, 
scalar quantities. Thus 

(115) 0/ = Ui^ + Vi^ + Wx^ C/ = U/ + V/ + W/. 

The mean energy of peculiar motion per molecule is ^mG\ and we shall write 

(1-16) ™.C? ^ I = 3BT., „.C? ^ I ^ 3EI., 
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where R is the '' universal gas constant/' These equations define /i and T ; we shall 
call Ti, T2 the absolute temperatures of the component gases. According to the 
theorem of equipartition of energy, in the uniform steady state of a gas they are 
equal. In the slightly disturbed states which we shall consider, the differences 
hi—h^, T]— T2 will be small. '^' 

We further define p^^ p2^^ Poy the mean hydrostatic pressures of the separate 
components of the gas, and of the total gas, by the equations 



(1*17) p, = ii^imjC/ ^ ^ "^ R^iTi, P2 = ^^2^^202" ^ ^ "" ^i^2^2. 

(1-18) ^.^p,+^. = i + ^^^A^B,T.. 

The last equation also defines h^ and Tq ; the latter will be called the absolute 
temperature of the composite gas. Clearly 

(1181) To = A,T, + A,T„ ~-=^ + ^- 

ho h^ k^ 

We shall define T'o^p'o? ^^0 by the equations 

(1-182) T'o = A,(l\-To) = -X2(T.-To) = XA2(T:-T,), 

h'^ ^\hi hj \h^ hj ^ '''Xhi h^t 

(1-184) y„^_|_^ = XA.(g-^^) = E.oT'„. 

The following equations are inverse to the above : — 
(ri85) T, - To+^% T, = To-|-^ 

Ai A2 

(1*186) pi ^ X2{Po-^hp'o)^ P2 = ^i(Po~Wo)- 

^ [In the paper as originally communicated, no account was taken of these differences, a preliminary 
examination having indicated, that they do not materially affect the theory of diffusion. The distinction 
between Ti and T2 has been re-introduced at the suggestion of a referee, in order that its influence, if any, 
on the phenomenon of thermal diffusion might be made clear. It will appear that Ti - T2 is a small 
multiple of the rate of change of A'o (or Xi) with time, so that in steady states of the gas it is a small 
quantity of the second order only 3 in particular, the phenomenon dealt with in § 16 is unaffected. 

Throughout the paper, wherever a distinction is made between hi and Ag, Ti and T2, or wherever T'o 
(equation I '182) appears, this has been introduced on revision (June, 1916). The original form of the 
equations may be found by making the difference zero. An appendix has also been added on acconnt of 
this extension. — June ^, 1916.] 
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(d) Some Convenient Limitations of the Problem. 

Our immediate aim is to determine functions /i(Ui, Vi, Wi), f2{V2, Vg, Wg) which 
define the distribution of the peculiar velocities of the molecules, i.e,^ which are such 
that the number of molecules of the group considered (the appropriate suffix 1 or 2 
being added throughout), the components of whose velocities lie between (U, V, W) 
and (U+<iU,V4-(^V,W + dSN) respectively, is 



(119) 



vfiU. V, ^N)dUdyd\N 



per unit volume. Besides the independent variables U, V, W, these functions 
/(U, V, W) will also involve such quantities as v, X, P, c, h and their derivatives, all 
of which are functions of (^, y, z, t). The distribution of the peculiar velocities is, 
however, clearly unaffected by the absolute magnitude of the mean velocity Cq 
(though the same is not true of the derivatives of Cq), We may, therefore, legiti- 
mately make the simplifying convention that c^ = at the particular point and time 
under consideration. This merely amounts to a particular choice of uniformly moving 
axes of reference, a choice which the laws of dynamics leave quite unrestricted. 

Our concern being with problems of molecular rather than mass motion, we shall 
suppose that the acceleration of the gas as a whole is of the first order only, which 
requires that the resultant force on unit mass of the gas, viz., Pq/w^o? shall be small 
We shall also suppose throughout that the velocity c^q of interdiffusion, and the 
derivatives, with respect to space and time, of i/, X, c, h, are all of the first order, at 
most"^ ; consequently, since in this paper we shall neglect second order quantities, 
products and derivatives of any of the small quantities just mentioned will be omitted 
from our analysis. 



§2. The Equation of Teafsfee op Molbculae Peopeeties. 

{a) The Equation of Continuity, 

The general equation of transfer for a function Qi of the velocity components 
{u)i, {v)i, {w)i of a molecule of the first kind ist 



(2-01) 



AQi = TT (flQl) + s 



dt 



z,y,2 



't'^^^'i-^Mw\. 



where AQi denotes the rate of change of j^iQi at {x, y, z, t) produced by the encounter 
of the molecules of the first kind with others of the same or the other kind. 

■^ [And likewise /oi T'o, l/^'o, when we are considering unsteady states in which Ti ^ T^.'—Jwm S, 

J. if JL vj • I 

t Gf, Chapter IX. of Jeans' 'Treatise' (2nd ed.), and also, for the details of the reduction of (2-01)- 
(2-02) and (2*09), ' Phil Trans./ A, vol 216, p. 285. 
VOL. CCXYII. — ^^A. T 
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If we write Q = 1, since AQi = 0, and . /^\ = 0, the equation of transfer reduces 



to 
(2-011) 



dt dx dy dz 



which is the equation of continuity for the first component of the composite gas. By 
adding to this the corresponding equation for the second component, we obtain the 
equation of continuity for the gas as a whole, in the form 



3 (^1 + ^2) , '^ {viUi-\- V2U2) , 3 (i^i^\ 4-1^2^2) , 8(1^1^14-1/2^2) _ 



dt 



+ 



dx 






dy 



+ 



dz 







or 



(2*012) 



dt dx dy dz 



0. 



If at the point under consideration the mean velocity of the gas is zero, the last 
equation may be written 



(2-013) 



1 ^ 

Vq dt 



dU() , ^^'o , dw^X 



+ TT-^H- 



dx dy dz I' 



(b) Qi = UiCi^ 

If in the equation of transfer (2-01) we assign to Qi the value {u)i{c)i^% and omit 
all terms which in a gas of ordinary density are of the second order, we find that 



(2*02) AUA'' = 1.3.5...(2g + 3) 



i^i 



{2him,y 



+ 



1 3 



i^i 



dUQ 

. 3^ ' 2himi vi dx nil 



n. ox\2fiimi. 



We have here used the convention that Cq = at this particular point and time, 
so that (c)i and Ci are identical, and, except in difierential coefficients, Cq can be 
neglected. 

If we multiply both sides of (2 '02) by mi, and add the corresponding equation for 
the second group of molecules, in the case when 5 = we obtain the result 



(2-03) A (miUi + maUs) = (vimi + v^m^) -r-^ + 

ct 



duQ , / 1 3i/i 



+ 



1 3 



^2 



2hi dx 2^2 dx, 



- (.A + .,X,) + |''i^(^) +". ^(^ 



dUr 



3 / 



/ ^1 



^2 



= .0^0 -^» -.0X0+ 3^^-;+ 2^ 



= v^rriQ 



dt 
du 



dt 



^ — ^0X^0 + 



3po 

dx ' 
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by (118). In this equation the left-hand side is the rate of change of the momentum 
per unit volume of the whole gas at {x, y, z, t) which is produced by encounters 
between the molecules. Since, however, an encounter between two molecules leaves 
their combined momentum unchanged, A (wiUi + mgUa) is zero, and consequently 

(2-04) ^ = J_X,— i-l^o. 

ot tUq vqITIq ox 

This and the two similar equations in y and z are the equations of mean motion of 
the ffas. 

o 

We will now apply (2"04) to the elimination of -^ from (2*02). At the same time 

we shall neglect the difference between hi and A-o? Tj and To in products or 
derivatives, since our equations are to be carried only to the first order of accuracy. 
Then (2*02) becomes 

(2-05) 3 i2h,miy ^i ^y g 2. ^ 1 Idv^ ^m^^ ^ ri!±^^(^^i)±(l. 

1.3.5...(2§ + 3) 1/1 ^ ^ 2^0*^1 3^ ^0 ^ ^ v^niQdx dx\2hQ/ 

2hQ\vidx vq ax/ dx\2hQ/ Xi v^^mQ dx 8x\2% 

aT. 



= L 

Thus, if we write 



1 3X^0 ___y/ _ m^Q 3p(i 
2^0 3^ ^ ^o'^o 3^ 



+ E5 



dx 



(2-06) ^'^= 1 ^_^, m^o^^^ 

2ho dx vqWIo ax 

the equation (2*05) and the corresponding equation for the second set of molecules 
become 



To this order of approximation, therefore, these equations do not involve TV 

(c) Qi = Ux^C,^ 

When Qi = UiGi^% the equation of transfer takes the form (c/. § 3 (C) of my second 
paper, loo. cit.) 

(2-09) 15(2A,m,^^ .^J.c■>^ ^ 5 Jlh.+ {s + l)2h,^(^) 

^ ' 1 .3.5...(2« + 3) ' ' 'U 3« ^ '' dt\2hj. 

\ dx dy oz / 
T 2 
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Now 
(210) 



O r t- - I — - — 

OX oy dz 






, '^ ^0 , '^ ^^0 \ I 2. / 2 A 



dVn dw, 



'^^0 



dy dz 



5 k^ jL^r 
Vq ut 



by (2*013), where also we have used the notation indicated by 



(2M1) 



^xx ~ 


_ ^ Buq dvQ dwQ 
dx dy dz 


^yy 


__ p dvQ dWo dUo 
dy dz dx ' 




__ ^ dWo BUq dvQ 
dz dx dy ' 



Evidently we have 



(2-12) 



Ca;j; + Cyj, + Cj,^ "• 



By substitution from (2"10) into (2 '09), we find that 



(2-13) 



15(2y^omi)'+^ 
1.3.5 ...(2s + 3) 



AUi'C/' = 5vJ- ^ + {s+l)2h 



(."1 



dt 



' dt \2hJ 



+ f (2s + 5)t/iC.^^. 



i(2. + 5)l| 

Vq Ot 







If we divide both sides of this equation by 10 . 2hi, and assign the zero value to s, 
it becomes 



(214) 






On adding to this the corresponding y and z equations, the result is 



(2-15) 



Ajm,0/ = fi>x^^-^^^-^'^-^ 



\Pi dt 



% 



dt 



There is a similar equation for the second set of molecules, and by addition of the 
two we get 



(216) 



A(i™,C.'+J»A') = = 4p.(i|-.-j L|), 



the left-hand side being zero, since energy is conserved throughout molecular 
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encounters. The integration of the last equation gives us the equation of state of 
the gas, viz., 

{2*17) i>o^o~''* = constant, 

or 

(218) T,vo-''' = constant. 

We will now return to the general case of (2*13) ; on neglecting the difference 
between hi and h^^ T^ and To, in derivatives and products, the equation becomes 









'XX 



by (2*18). The similar equation for the second set of molecules is 

(^■^«) 1 . \ ^1'%^ 3) ^"-^^-^ = -° IT ^^- ^'^"-'^ ----- 

From (2'20), by transformation of rectangular co-ordinates, or by direct calculation, 
the equations of transfer in this case may readily be shown to have the form 

(2-22) __|i|^!ig^AV,W.C,- = A(2s+5).;^,„ 



where 






(e) Q = C^'. 



By the addition of the x, y, and z equations corresponding to (2'19), (2'20), and by 
changing s to s— 1, we may obtain the following simple equations : — 

^^^"^^ 1.3.5...(2s + l)^^^ ~''~dt -~ 1.3.5. ..{2s + l)^^' • 
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§ 3. The Velocity-Distribution Function. 

(a) The Form of the Function. 

In the uniform undisturbed state of a gas, in which Cq, X, h are constant, while c'o, 
h^Q, P are zero, the velocity-distribution functions assume Maxwell's well-known 
form 



(3-01) 



(/i)o 



^0^1 \ '^-Aom,C,2 /y-X _ /^O^Y'^-Aom^Ci^ 



TT 



X 



The suffix appended to fi and /^ is to indicate the special state to which these 
equations refer. They clearly satisfy the necessary conditions 



00 



CO 



(3-02) 



(/)o c^U, dW, c?Wi = 1, j^^ {f,)oG,' dU, dV, dSN^ = 



00 



~ 00 



2^0^! 



and similar equations with suffix 2. 

In the general slightly disturbed state considered in this paper, /(U, V, W) will 
differ from {f)^ by an amount of the first order. From the equations of transfer (§ 2) 
it may be deduced^ that / may be expressed as follows : — 



(3-03) /i(U„V„WO = (/i)c 






Bx 



dy 



dz 



-^\2h,m,G, (c,,U,^ + c,,V/+c,,Wi='+2c,,ViWi + 2c,,WiUi + 2c,,UiV0 H^ (C/) 



-D, 



3t 



Ji (Ci^) 



(3-04) /^(U,, V,, W2) = U), 



1 -i2h,m,A, ( U,f 'o + y,^,', + Wa^o) F, (C/) 



-i2Km,B, (u, 1^ + V, |^° + W, ^"j G, (C/) 






The constants Ao, Bo, Co, Do, and the functions F(C^), G(C^), H(C^), J(C') 
remain to be determined. The latter involve {x, y, z, t) only through the occurrence 



The argument is given in my second paper, loc. cii., ^ 2, 6, and will not be repeated here. 
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of V and h (but not their derivatives) as factors. We suppose them to be capable 
of expansion in power series as follows^ (see Note C, p. 196) : — 



00 



(3-05) Fi(Ci^) = 2 a 



{2hmiy 



CO 



'■l.3.5...{2r+3) 



0,^, F,(C/) = 2 a_, 



{'ihm^' 



r = 



1.3.5...(2r+3) 



^^2 ? 



00 



(3-06) Gi(C,^) = 2' A 



(2Ami)' 



00 



r = 



1.3.5...(2r+3)r 



Ci^-, G,(C/) = 2' /3. 



(2/1^2)"" 



r = 



] .3.5...(2r + 3)r 



n 2r 



2 ? 



00 



{2hmi) 



00 



{2hm,^' 



(3-07) H,(C.)=E,,^^^=||^C,>, H.(C..)=2^,.,^^^^^^^^^|^C 



2r 

2 ? 



00 



(3-071) Ji(Cx') = 2 ^, 



(2Ami)' 



00 



''1.3.5...(2r + l) 



Cl ''5 J2(C!2 ) — S 0. 



{2hm2y 



r = 



1.3.5...(2r+l) 



p 2r 



^3 



The dash (') after the sign of summation in (3*06) is used to signify that the 
factor r in the denominator of the numerical coefficient is to be omitted in the first 
term (r = O). The choice of the notation +r and —r for the suffixes has a convenience 
which will become apparent later ; we may remark, in passing, that for this purpose 
a distinction must be maintained between + and — 0. 

In (3*03) and (3*04) the constants Aq, Bo, Cq, Do can be chosen arbitrarily, but 
when this has been done, the remaining constants a, ^, y, S all become perfectly 
definite. 

(b) Relations between the Coefficients. 

The velocity-distribution functions /(U, V, W) must satisfy the three conditions 
expressed by the equations (with appropriate suffixes 1 or 2 throughout) : — 



(3*08) 



f{u,y,\N)dudydys/= i, 



J J 



(3-09) 



jj[/(U, V, W)C'(^Uc?Vc?W = 



3 



2hm 



{hi or Ag), 



(310) 



/• /• 



/i (Ui, Vi, Wi) Ui t^Ui dVi cZWi = Ui-Ua = w'o/Xi, 



j|j/.(U. V. W.) U.<^U.c^V.c^W. = u,-u, = -u'JX,. 



* [Here, and throughout the remainder of the paper, where h and T appear without any suffix, they are 
to be read as ho and Tq.— /twie 2, 1916,] 
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These yield the equations 
(3-101) S^, = 0, 







oo 



(3-102) 2<5_, = 0. 







{3-103) iJ), ^ S (2r + 3) 4 = - 4f ' 

(3104) |Do%^S(2r+3)a_, = -ij^. 

Ot Agio 

(3-11) -i Ao^oSa^+Bof^SV-^^i = </A„ 

I do? J 

(3-12) -i Aof'o2«_, + Bo|isV-"^_i=-ttVx,. 

L da? J 

In each of the pairs (3*103)-(312) we will multiply the first equation by Xi and the 
second by Xg respectively, and add. We may separately equate the coefficients of 

^^Q and — (in the second resulting equation) to zero, since these quantities are quite 
independent of one another, and of their coefficients. We thus obtain the conditions 



00 • . 00 



(3-121) \i-E{2r + 3)S, = -S'o = ~\2^{2r+3)§ 



or, by (3-101), (3-102), 



(3-122) 2\iXrS, = -S'o = -2\,'SrS 



00 00 

r' 



(3*13) Xi S ex^ = — exo = — X2 2^ 



1 1 

00 00 

/v 

r? 



(3-14) Xi X'r-% = -l3'o = -As 2V-i/3 





00 00 







where also we have introduced a convenient notation for the separate sums involved. 
Expressed in terms of this {3'10l)~(3*12) are equivalent to 

(315) U% = l-fa^oAof ^o+^'oBo ^ , 

ot 

[Throughout the remainder of the paper we shall neglect -y^, ie., we shall practically 
assume that the ratio of mixture is not varying with respect to time. The values 
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of a, /3, y determined In the succeeding sections will, howeyer, hold good also in the 
general case ; the determination of the ^'s will be found in the Appendix. . All the 
corrections to (/)o are of the first order, and hence are separately deducible ; 
the complete value is obtained by adding to the series in a, /3, y, calculated in the 
body of the paper, the series in §, determined in the Appendix. — June S, 191 6 J] 

§4. Completion of the Equations op Transfer. 

(a) The Values of AQ. 

We may now complete the equations of transfer (2*07)--(210) by the insertion of 
the values of AQ, calculated in terms of the constants in the expressions for /(U, V, W) 
given in § 3. The calculation of AQ is a lengthy and elaborate operation which will 
not be described here, since a full account of it is to be found in my second memoir 
('Phil. Trans.,' A, vol. 216, §7, p. 301). It appears that only those terms in /which 
are of odd degree in U, V, W contribute to the resulting expression for AUC^*, and 
only the even terms, similarly, contribute to AU^O^*. The following results will be 
quoted forthwith : — 



(4-01) ^ilifef-miAUxC- = i 



mi ( N„ Aof Vr + sN'„Bo -^ M {pn (n«i) + pu (nsi)} 



(4-02) ^Jig^l-^,AU,0/' = ^ i [mi (N,Aro«. + ^N'„B„|f /3.)p.i (ns,) 

+ f^i \^„ Kt'^^-T + sN'^jBo g^ /3_ J {^22 (52*2) + /'2] {r^s^ } 



45 (2Am,)'+^ 1 A,,2n2«_ 1 



CO 



(4-04) ;:,\T':t ^AU/C/' = iC,c.. 2 N';,[yy,(rA) + y_.{Ate)+A(r2.2)}]. 

(b) Explanation of the Notation. 

In the above equations r and 5 may take all positive integral values, including 
zero. We shall presently alter our notation so as to consider also negative integral 
values, but in the following definitions of N^^, N'^^, N^''^^, the positive numerical 
values of r and s are in every case to be used on the right-hand side : 

(4-05) N„ = {2'•-»+H^+fU:(s+lX.l}-^ N',, =;;iN„, N"„ = N,,i„,, 

In the ease of N'^,, the factor r is to be omitted from the denominator when r = ; 
when 5 = 0, the value of sN'^„ in which form N'^, occurs in (4'0l), (4*02), is to be taken 
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as r~^N„. The meaning of the symbols (r+f)^^.!, (s+f),+, will be understood from 
the following definition, where q is a positive integer and p any number whatsoever : 



(4-06) 



i? =^(p-l)(i>-2)...(_p-g+l), Po = l- 



We shall have frequent occasion to use this factorial symbol. 

The following are the expressions found for p (rs) and p' (rs) in the various 

C9jS6S I ■"" 



(4-07) 



pn (riSi) = ^hi 



g-fa'+y')a,2^2 






+ 2y'm^ B'^^'ir, s)+-^B''-'{r, s)-B''{r, s) 



(4-08) pi^{riSi) = 



■hffiVg 



14^+1 



4^+1 



dxdy. 



(r+l, s+1) 
il; = 

+ AfjL^f^yjr^^^ (y) B" (r, s)']rm dx dy. 



(4-09) pi2 {r^i) = ^vii/gMi/ 



I, 



g-fa^+j,»)^2^2 



(f + 1, S + 1) 



2 ( - 1 )* [<p,J' (y) {/«,;'** (r + 1 , s) - 2 (M1M2)- VB* {r, s) 



k = l 



+ M23*(n s+ 1)} -4 (/*i/*2)V^n2*(?/) B^(r, s)l„2 dxdy. 



(410) />'n (nsi) = ¥"1' jj 



g-(x^+y')a;y 



4(r+2, s+2) r 

^' <f>u''{y) B'''{r+2,s) + Wir+l,s+l) + B''{r,s + 2) 
1 L 

+ 4?/^ |||±y('b^*+ H*- + 1 , «) + B'*+ Hn s + 1 ) 

B^*(r+l,s) + B^*(r,s + l) 



+ 



2k 



ik + 



JB'''-'{r+l,s) + B'''-'{r,s + l))\ 



f(2^+2)(2^ + l)g.,,.( ^^ 



+ 



.(4^ + 3)(47£+l) 
(2A;+1) 



{2hy 






4. 2a:(2a:— 1) -psk-i (^ ^\ 
^{4k+l){ik-l) 



2,'|±1b-h.,.) + JLb-(^,^) 



dxdy. 
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(4-11) /.^(nsi) = !-.,. J[e-<^'«'>a^2/=' 



•/ «/ 



ir+2, s+2) 

2 li>i2Hy) {B*(r+2, s)+|B*(r+ 1, s+ l)+B'ir,s+2)} 

+ Bfii^y^^ii {y) {B^ (r + 15 5) +B^ (r, s+ 1)} 

•.t 

- 1 &fxiy%i{y) B* (**, s)]i2i2 dx dy. 



(4-12) ^(^a) = i«'i''2 e-'^^+^^'^y 



V V 



Cr+2, «+2) 

+ 4Mi/i2f 'B^ {r, s)}-8{fx,fA,yf^firJ{y) {fxJ^^B'ir^ 1, §) 

^2{^,M,)"VBMn^)+/x,3Mn^+i)} 

There are also six other equations^ similar to the aboYe^ except that the suffixes 1 and 
2 are interchanged ; these need not be written down here. 

The limits of integration of x and y throughout the above expressions are andoo. 
The upper limits of the summations are in each case indicated by two numbers, which 
are not necessarily integers ; the upper limit is to be taken equal to the greatest 
integer which does not exceed either of these two numbers. The suffixes 1212 or 
2112 on the right-hand of the main square brackets of (4*08)^ (4*09)5 (4*11), (412) are 
there placed only for conyenience in printing : they should really be appended to each 
of the symbols B^' {m, n) contained within the brackets. These symbols are defined 
by the following equations : — • 

(413) B^W2 (^. n) - ^A'i2f,,x% 2m#) . **A*(2pt#, 2ja,y% 

where '^A!'{u^ v) is a polynomial in powers of u, v defined thus :— 



■a' 

(4-14) »A" {u, v) = ™A* (v, u) = (-Y S jf\. ^?^f iv' «'""''^* k^m. 

When h > m,'''M {u,v) is zero. Also mb M2. ^12? ^21, haye the following values :— 

(415) Ml = mJimi+ms), fi^ = m2/{mi+m2% M12 = mi^2 = Mi//^2, /^2i = W^i = /"^/mii 
so that 

(416) Mi + iW2 = 1? Ml2/^21 == 1' 

u 2 
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Finally we must define the functions ([J' {y), xj/ (y), x^ (p) I these are the only 
quantities in the expressions for p (r, s) and p (r, s) which depend on the law of inter- 
action between molecules. If two molecules mi, mg encounter one another, the 
direction of their relative velocity will he changed through a certain angle $12 in a 
plane parallel to the initial relative velocity and to the perpendicular (of length p) 
between their initial and final lines of undisturbed rectilinear motion ; denoting the 

magnitude of the initial or final relative velocity by ( -~ ^ y^'y, 0^2 will be a function 

of p and y, the nature of the function being determined by the mode of inter-action 
between molecules in proximity to one another. Then we have 



.on 



(417) ^1/ {y) = {2k + 1 ) {mW^ (wi + m,) } -''-y { 1 - P^ (cos ^i^) } p dp, 



.00 



(4-18) i^vHy) = {2k + l){fA,fA2h{'m,+m2)}-'^j (l-cos 0i2)Pj(cos Oi^jpdp, 



.00 



(4-19) xi3*(?/) = {2k+l){fx,fx.Ji{m, + m,)y^-y (l-cos 6,,yF,{G08 e,,)pdp, 

Jo 

* 

where P^. (cos Qi^ denotes, as usual, the Legendre function of cos 6^2 of order h. By 
changing the sufiix 1 or 2 throughout into 2 or 1 respectively, we obtain the 
corresponding expressions for ^^i or (p^i, and so on ; mere interchange of the sufiixes 
does not aifect the functions. 

By means of the recurrence formula for the Legendre functions, viz., 

(4-20) (^+l)P,^i^(2^4-l) cos 0i,P, + /cP,„, = 0, 

we may express yfr and x in terms of the function ^. Thus for yf/ (y) we have 

(4-21) riy) = j^i>''Hy)-'PHy)+^^<p''-Hy). 

In this way we may prove that 

(4-22) i.^y) = i^^y\ ^ (y) ^ §i>^ ^-^,i^ {y\ ^ iv) == ' <pHy) + UHy)^ 

From the symmetry, with respect to r and s, of the expressions on the right hand 
of equations (4-07), (4*08), (4*10), (411) it is clear that 

(4-23) pn(n^i) = PniSinl Pn{nSi) = pn{s,n), ^12(^^1) = ^12(^1^), /i2(n^i) = pAsiri). 

(c) Special Values of p (r, s) and p (r, s). 

To facilitate the exposition of subsequent parts of the work it is convenient at this 
stage to write down certain special cases of the equations (4*07)--(412), after executing 
the integrations with respect to x and y. Owing to the generality of the functions 
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(j), Yy, X (depending on the functional relation between 0^2 and p, y) the integration 
with respect to y can only be made formally, however, and for this purpose we shall 
use the notation 

(4-24) KMO^/,A",x e-^y(?/)?/^^*"'-^>^2/. 

The suffix 11, 12, or 22 is to be appended to K^(^) to correspond with the suffix 
of 0^ {y) on the right. The notation is chosen so as to make K^' {t) equal to 0^ {y) 
when this is independent of y^ as in the case of molecules which obey Maxwell's 
fifth-power law."^ 

For small values of r and s it is convenient to simplify our formulae by writing also 

In terms of the above notation we may now give the following results : — 



s-l 

S 





4-26) pn (nOi) = /5h(0a) = 0, />n(Sil,) = Pn{liS,) = l%irvMs + il^i^'M ^s-iGtk\„ 

4-27) p,, (siOi) = p,, (OiSi) = i7rv,v,fji,2' {s + |),^i K'i2 (0) 2, G,^rWk, 



4-28) p,, (r^Oi) = - |7r.i.2Mi2'- (r+f ),^i K\, (O) 2rC,M/-Vi'^t, 



4-29) p,, (O^Si) = — |7r.i.,Mi2' (s+f X+i K'j2 (0) 4c,Mi'-VA, 







4*30) pi2{Sili) = PuihSi) = |7n/ii/2M22'(5 + f ),+iK^2(0) f 2,C,m/"W U^^ + Ma^^+i) 



L 



s-l -| 



4'31) pUUSi) = — |7n/ii/2Mi2^(5 + fLiK'i2(0) |2,C,M/-W(M2^i + MA+i) 





s-l -| 

+ ^/^lM2 -^s-l^itMl M2 V^^« + ^<+2"~"5^32 > 

J 

4-32) /n (Oa) = /n (^lOi) = ^5 W (^+f X-^2 K\, (0) 2,CA„ 



4-33) /i2 (Oi«i) = p'u (sA) = f i7r.i.2M22» («+D.+2 K'l, (0) i;,o,M/-v/ (mi^*+ Am^^W. 



4-34) p\, (0,Si) = f|7r.i.,M22' (s + |),+2 K'i2 (0) 2,Qmi'-V/ (-Mi^* + tVi^*12). 



* (7/ * Phil Trans.,' A, vol. 216, § 9 (0), p. 323. 
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It may readily be deduced from the above equations that 

(4*35) pi2 (nOi) = -P21 (nOs), P12 (^20i) = -P21 (nOs), 

(4*36) /^iPi2 (Oa) = — /^2Pi2 {O2S1), P1P21 {O1S2) = — /a2/)2i (02^2). 

Also, in the case of Maxwellian molecules, for which k^ = 1 and ¥12 = ^^ 
(whatever the value of t) we shall have 

(4-37) pi2 (ll5l)-5/)l2 (Oi^i) = %7rP,P2f^2^'S (5 + -|),+i K^32 (O) (3/// + /X2' + 5MlM2^'l2), 

(4-38) /)i2 {l2Si)-5p,2 {02S^) = - i7r.,.2MiM2'2^^ (^ + t),^i K'12 (O) (4+ P^. 

§ 5. The Symbolic Solution por the Coeeeicients in / (U, V, W). 

(a) The Linear Equatioiis for a^ /3, y. 

We now refer back to the two corresponding sets of equations (2*07)-(2*10) and 
(4*0l)-(4'04). For the two members of each pair of corresponding equations the left- 
hand side is the same, so that we may equate the right-hand sides. Also, as regards 
the first two pairs, we may separately equate the parts which contain as factors the 

independent quantities ^\ and y- . Thus we have 

-i 00 
(5*01) - Ao 2 N,,[mia^{/)a(n5i) + /)i2(n5i)}+^2a_,pi2(r25i)] =1, (5 = to 5 = 00), 

1 00 

(5-02) Ao 2 N„[mia,^3i(riA-2) + OTa«_,{/)22(^2*'2)+/>2ife)}] =1, (s = to s = oo), 

I'd r = 

(5-03) :5-Bo S N',,[m,/3,{/)n{n5l)+Pl2(n5l)}+^2iS_r/>12(^'25l)] = 1, (5 = 1 to 5 === CX)), 

lSiV\ r=:0 

~i 00 

{5'04) :|5~ Bo 2 ^\^Xj^Srp2l{n8^-^m2^-v{p22{r^^^^^ =1, (5 = 1 to S = CX)), 

■(5-05) ^s- B, 2 sN'^,[mi/3,/5j2(7-iSi)+m2/3_,/5i2{^2Si)] = 0, (s = O), 

1 00 

{5*06) -Co 2 N'',,[y,{/ii(n'^i)+/i2(n5i)} +7-^/12(^^1)}] ==1, (5 = to 5 = 00), 

(5-07) - Co S W\, \y,p\^ {r^s^-^y-r {922 {^2^2)^ p\i (^^2)}] = 1, (^ = to ^ = 00). 

By virtue of (4*26) and (4*35) the equations (5*01) and (5*02) are identical when 
5 = 0, while (5'03) and (5*04) then assume the same special form (5*05) ; we may 
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recall that in § 4 (6) the value of sN'^^ ^^^ defined to be r""^N^g when 5 = (or N,,^ 
when r is also zero). Moreover, from (4*26) and (4*36), it is clear that in (5*01)-~(5*05) 
the ratio of the factors of cto and a_o, or of /5o and ^_o has the constant value — 1. 
Consequently, the above equations do not enable us to determine the separate values 
of ao and a^o? ^^ A and /5_o, but only of aQ—a_o and ^o~"^--oj which form single 
unknowns. When these and all the other values of a^ and /3^ (r 52^ O) have been 
determined, the separate coefficients cto, a__Q, /3o, /3_o may be deduced with the aid 
of(313)and(314). 

In order to simplify the notation of our formal solutions for the coefficients a, ^, y, 
it is convenient to re-write (5*Ol)--(5*07') in the form 

(5*08) S a^^a^+aQs{aQ—a_^Q)+ S a^,a^ = 1 (^=— ooto5=oo, including 5 = O), 

r= —CO r = 1 

~l o^ 00 

(5-09) 2 6r+i,»+]/3r-&»TT' + ^ K-i,.-A = 1 {s= -co to s = 00, excluding s = O), 

-1 Of ^ 

(5-10) 2 &r+A-&^+ S &r-i/3. = 0, 

r = — CO ^1^2 T — 1 



CO 



(511) S Cy^yr = 1 (±^? ±^ range from to 00), 



y = — 00 



In obtaining (5*09), (510) from (5'03)-(5'05) we have eliminated ^^ and ^_o ^J 
means of (314), and in order to preserve symmetry we have subtracted l/^Xi or 
1/5X2 times the equation (5*05) from (5*03) or (5*04) respectively. The new symbols 
are defined as follows : — 

(5*12) 5>0 ao, = -a^Q, = — AoNo,/)jt2(Oi5i), a^, ='—^^tApn{nsi)+pn{riSi)) r>0, 

(513) S<0 ao,=~a„o.=~— AoNo,/>3i(Oi52), %, = ~^AoN^,{p22(^2^2)+P2^^ ^<0, 

(5*14) r <0, 5>0 a^, = ^AoN^.p^^i'^^^,), a,, = - — AoN^^psi (n^s) ^^> 0, 5<0, 

(5*15) r^O b^ =■ -^ (^r+l.O^^OO/j ^ = A -o <^00? ^r = A -o' , I (^r--l,0"~^--00/ ^^— 0> 
(olu) 5]^0 0= — ■ •— — — -■ ^ TX^otJLi — ^,10/, Q z=i -^ — ■— — ^r ""~~~r 1^0 -_i — ^00/ 5*^ — 0, 






(5-17) s>0 ^ 



1 T> 1 

^" " Xl A^R (r+1) (s + 1) ^«^'-M+i-«-o,.+i-«r-i.o+«_oo} ^<-0, 
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{5-18) s<-0 -( 



^'■'~ ~X,A^{r + l){s + l)^'''^'"-' '"'"-' ^-+'-o + ^««^ *'>^' 



Ks = - 



1 B 



A. 



a;;b;{^+i)(5+i){^'-i--i ^-"'-i «.-i.o+«-oo} ^<-o 



(519) r>0, s>0 c„ = -GoW',, {p'n (nsi) +/i2 (nsi)} 



c„ = - CoN"„p',2 (r^Si) r < - 0, s> 0, 



(5-20) r< -0, 5< -0 c„ = i CoN",, {p',, {r,s,) +p',, (r^s,)} 

Cr, = - Go'N"rsp2A (n«2) r > 0, S < - 0. 

On the right-hand sides of the above equations the positive numerical values of 
r and s are to be used (whatever their signs on the left-hand side) except when they 
are suffixed to a or 6. 



(b) The Formal Solution of the Equations for a, ^, y. 

If we may solve the linear equations (5'08)-(5'll), each containing an infinite 
number of unknowns, as if they were finite, we arrive at the results 



* \affin) 



oto— a_o 



a^ = 






(r = — oo to r = ooj excluding r = O), 



A 



= ''/) "^"[^ (^r = — CO to r = ooj excluding r = O), 



^ \^mn/ 



K 



^' V (c„„) 



(-7' = — oo to r = + oo J including r = ± O). 



mnf 



In these equations V (a^„), V (&^„), V (c^„) denote the determinants which have a 
^mny ^mn ^^ their general element ; in the two latter, ±m^ ±n range from to oo^ there 
being also a central row and a central column in V (6^^) which are not enumerated 
by m or n. In V (a^„) the values ± of m and n are not distinct from one another, 
so that this also has a central column {m = O) and central row (^ = O) ; V (c^„) has 
not got either of these, since ±0 correspond to different rows or columns. The 
determinant V^ denotes that obtained from the corresponding V by replacing all the 
elements of the r*^ column by unity or, in the case only of the central element of the 
r^^ column of Vy(&^„), by zero. It may be remarked that all these determinants V 
and V^ are infinite in both directions, covering the whole plane. In two quadrants 
(m, n both positive or both negative) the determinants V possess symmetry 
(cf (4-23)). 



MON ATOMIC GAS: DIFFUSION, VISCOSITY, AND THERMAL CONDUCTION. 139 

In order to make this investigation complete, from a purely mathematical 
standpoint, it would clearly be necessary to supplement the above formal solution by 
a discussion of the questions of convergence raised in the course of our analysis. The 
complexity of the problem, however, and the rudimentary condition of the theory of 
infinite determinants, may well absolve the author from an attempt at such a task, 
for the present. From the standpoint of mathematical physics there is, fortunately, 
sufficient numerical, evidence {cf, § 136, jT, and ' Phil. Trans.,' A, vol. 216, p. 330, 
Table III. ) to afibrd reasonable assurance that our expressions converge satisfactorily ; 
this is especially so in regard to the formulae actually used in the applications of 
/(U, V, W), U., (5*22)-(5-25). 

(c) On Certain Combinations of the Coefficients a, ^, y. 

For the purpose of the theory of diffusion, viscosity, and thermal conduction in 
composite gases, we require only certain combinations of a, ^, y, and never their 
individual values. The following expressions comprise all those we shall find 
necessary in this paper; in connection with them we may refer back to (3*13), (3*14), 
and the formula (3*03), (3*04) for /(U, V, W) :— 



c» 



(5*22) a'o = — XiXgS (a^— a_^), 







c» 



(5-23) /3'o = -XiXJ(i8o-j8_o)+2r-H^,-^_,) , 



c» 



(5-24) 2;i,(..imi3Ui'-CiH..jm23U/-C/) = -^^Q,c^,^{viy, + v^y_X 







(5*25) 2h{vvmMS^x^-\-v,m^\Jj:^i) = -f Ao^'oi 1(^1^0 + ^2^.0) + S (r + f) (i/ia^ + i/^a.,) 







We therefore desire to obtain concise expressions for the following quantities : — 

CD C» 

(5-26) 2 («,-«_,), /3o-/8_o+ 2 r'' {^-^-X 

1 

c» 00 c» 

(5-27) i:r{v^a, + v^0L_,), ^{vSr+^^^-r), 2(^iy, + .'2y-r)- 

1 1 

The denominators in the expressions (5*21) for a^, ^^, y^ are independent of r. Hence 
our problem consists in the combi^a'"tioxl of the numerator determinants V^ for an 
infinite number of values of r, 
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(d) Difference- Transformations of Infinite Determinants, 

It is convenient at this stage to describe certain operations, by the application of 
w;iich we are enabled to preserve the symmetry and increase the convenience of our 
determinantal formnlse. These operations will be termed " continued differencing '' 
by rows, by columns, or by rows and columns, and we shall denote the corresponding 
symbolic operators by 4#j 4o? ^^ Ks respectively. We shall first define them in relation 
to infinite determinants which cover only a quadrant of the infinite plane (r ^ 0, 5 ^ O), 
and afterwards in relation to the more complicated type which occurs in this paper. 

The operation of continued differencing by rows {S^,) or by columns (S^q), applied 
to the determinant V (fs), where r ^ 0, 5^0, transforms it into the determinant 
V {4/„) or V (4„/„) respectively, where 

s r 

1^ ^"/ ^Osjrs = ■^ (—1; s^nJr,s~m ^rOjrs — ^ \ — 1 j T^mjr~-m,r 

We may effect the operation Sq^ as follows : from each element of row s we subtract 
the corresponding element of row 5 — 1, for every row from s — 1 onwards : this done, 
we repeat the operation on the transformed determinant, except that we now begin at 
5 = 2: and this process is continued without end, beginning each time with the row 
next after the initial row on the previous occasion. It may readily be seen that the 
result is as we have already stated, and that the value of the determinant is 
unaffected by the operation. Continued differencing by columns is strictly analogous, 
and need not be separately described. Continued differencing by rows and columns 
is performed by applying the two separate operations successively, the order being 
immaterial. Without alteration in value, the determinant is thus changed in form 
from ^ {frs) to V{4s/r^), where {cf (5 '28)) 



Iti o 



m ~ n = 



In the case of determinants which are infinite in both directions, the operation of 
continued differencing by rows and columns is effected b)^ applying the process 
described above to each separate quadrant ; thus the differencing by rows is performed 
by differencing outwards from the centre row in both directions (above and below), 
and likewise, by columns, both to right and left of the centre column. Neither of 
these partial operations, nor the complete process, alters the value of the determinant, 
a fact which we may express by the equation 

(5-30) V {/.) = V (5,„/„) = V {SoM = V (^„/„), 

where the notation is similar to that used in the' former case. There, however, r, 5, 
and consequently also m, n, were necessarily positive or zero, while our convenient 
notation for determinants of the present type involves also negative values of r and s. 
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In order to make the former definitions of S^^, S^^, §j.s applicable to this case, we must 
adopt certain conventions as to the interpretation of (5*28), (5*29), to allow for negative 
quantities. These conventions are {a) that in ^C^,sC„ the positive numerical values 
of r and s, m and^n^ are to be used in all cases : and (b) that r—m and s—n retain 
the same signs as r and s respectively even when m — r, n =^ s. The latter rule 
preserves, in relation to f^^^ the distinction between +0 and —0 which in certain 
cases we desire to maintain. 



00 



(e) A Symmetrical Expression for 2(a^-~a_^). 



In discussing S(a^— ot_^) it is convenient to change our notation for V^(a^^^)' 
writing it in the form V (^a^^) to denote the determinant whose (m, nY^ element is 
^a^^. If we define ^a^^ by the equations 



(5*31) ,a^^ zz a^n (^ ^ ^), r<^rn = Ij r^mO ~ r^m, -0 

we make V (^a^„) identical with V^(a^^) as defined in § 5 (6), and therefore, by (5*21), 

(5*32) ^^^^^^^pi^, ^^::=.pra^ ^^^^y 

By applying the operation of continued differencing by rows to V {^a,^^) we transform 
it into V {Son.r<^mn)^ whore 4n.r^m« ^^ defined by (5*28) (putting f^^ ~ ^a,^^, and making 
no distinction between /^o and /^, _o — cf. the third equation of (5*3l)). Now it is 
readily evident that 

(5*33) So^^^a^,^ ^ So^a^^ (m 9^ r), <^On.r<^r« =0 (^5^ O), (^oO.r^rO =^ 1, 

so that V(4n.r^mn)? ^ud consequently, also, V^(a^^), is equal to ^^r{^on.<^mn)^ this being 
defined as identical with ^ (Son.^mn) except that in the r^^ column the central element 
is unity, while all the others are zero. Hence '^^r{^on<^mn) is clearly equal to the minor 
of the r*^ element of the central row in V ((^on^mn)- If we replace the elements of this 
central row by +1 (m>0) or — l(m<0), and denote the result by V^ ((5o,,a^J, we 
may evidently write (c/. (5*32)) 

'* "'^ V (a^n) ^ (<^On^mn) 

by (5*30). 

We next apply to V {S^,,a.^n) ^nd V (So^a^^) the operation of continued differencing 
by columns, so as to transform them into ^ {§mn^mn) ^^^^ '^^ (^mn^mn)- We make one 
slight diff'erence of rule here, as compared with the former differencing by rows : that 
is, in (5*29) we shall preserve the distinction between m = and m = —0, writing/o^ 

* The signs of m, n are the same, of course, as those of r, s respectively. 

X 2 
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for the elements of the central column in V or V\ and defining /_on ^s equal to — ,/on- 
This means, in effect, that the elements of the central column are to be reversed in 
sign, for the purpose of differencing on the left (in < O).^ Now, by definition, 
^^ (^on^mn) differs from V (^on<^mn) ^^Y ^^ ^^^ central row, for which 

/^o = 1 (^ > 0), /^^ = - 1 (m < 0), and also /.oo ^ - 1, 
so that 

^mJmO = if m 5^ 0, (Joo/oo =1- 

Hence, when transformed into V^ ((^^^a^^)? it differs from V {S^^a^^ only in the central 
row, all the elements of which are zero except the central one, which is unity. In 
other words, V^ {^mn^mr) is the principal minor of V (<^^«a^), and the expression for 

00 

2(a^— a_^), viz., 



2 



(5-35) 2 («,-«_,) = Z?r''^i = - 6. 

has thus been reduced to a concise symmetrical form. 

{/) A Symmetrical Expression for ^V 
It appears from (5*21) that 

where Vq is the same as V, except that all the elements of the central column are 
unity save the central one, which is zero. If we transform V and V^ by the operations 
described in § 5 (d), differencing by columns with m = 1 (for the right) and m — — 1 
(for the left) as starting points, and similarly for the rows, we leave the central 
column (not enumerated by m) untouched as regards the first part of the operation, 
and the central row untouched by the second part. Thus we obtain the result 

(5-36) ^'o==-\y^2Z'l'^-i-}> 

V (4«6^ J 

where V(^^„6^^) is the determinant whose general element is SjnJ)^j,{±m, ±n ranging 
from to CO )^ and which has a central column Sorfin> ^ central row S„,ob^, and b as the 
central element. From this determinant we obtain V^^ (S^rfi^n) by substituting zero 
for all the elements of the central column except those next to the centre on either 
side (n = ± O) which are replaced by unity. 

^ We may note that the elements of the central column are the coefficients of ao in (5-01)-(5-05), while 
the elements reversed in sign are the coefficients of a«o- 
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OO 00 oo 



(g) Symmetrical Expressions for ^r{via^-^v20i^r)^ 2(i/i/3^4-i/2/3_r)j 2 (^ly^-fi/^y-r). 

1 10 

Ab in § 5 {e) we may prove that 

2.r{via^-\-v^a_^) - . v > 

where ^o{^on<^mn) is identical with V (5o„a^„) except in the central row, the elements 
of which are equal to rvi (r >> O) or rv2 {r ^ O), the positive numerical value of r being 
taken in both cases. On applying the operation of continued differencing by columns 
we obtain the result 



(5-37) ir{v,a, + v,a_,) = ^Si^, 

where Vo(^^„a^^) differs from V(^^„a^„) only in the central row, all the elements 
of which are zero save those on either side the centre {r = 1 and r ~ -~l), which are 
respectively equal to vi and P2' This follows from the fact that 

V = {r^ 1), V =1 {r^ 1). 

Again, from (5*21), we may prove in the usual way that 

where V^ and V are the same, except that in the r^^ column of V^ all the elements are 
zero save those on either side the centre (i.e., n ~ ± O) which are unity. If in 
V and V^ we add half of each of the rows 7^ = + to the centre row, and subtract 
this new centre row from the rows n — ± 0, V becomes transformed into V (^on^^'mn)? 
where b'^^ = 6^„, except when ?^ = ± 0, while 

Similarly, V^ becomes transformed into ^^((^on^^mn)? identical with ^ {^onb^mj)^ except 
that in the r^^ column all the elements are zero save the central one, which is unity. 
Consequently we may write 

2 1 Q A Q \ ^0 \^QrP mn) 

where Vq is the same as V except that in the central row all the elements on the right 
of the centre (m ^ O) are equal to j/j, all those to the left (m <^—0) are equal to P2, 
while the central element is zero. We now " difference ^^ by columns, with the result 

(5-39) SMr + ^^/S-,) = 54HH- 
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The general element of V is ^rnJ^rnn ^^ ^mrfi^mn^ which are the same except when 
^ — ± ; the general elements of these two rows are ^^o^^mo? <^m-o^^m-o (<?/• (5*38)), the 
general elements of the central row and column respectively are S^o^\^ and Sorfin^ while 
the central element is b. From this we obtain Vq if we replace all the elements 
of the central row by zero, except the two on either side the centre, which are 
p^ {m ~ O) and V2 {m — ~0). 

Finally, we may prove that 

[0 4:0) Z [v^yr'^P2y--r) - ry/. ^ V ' 

where V is defined as usual ( ± w, ±n ranging from to oo ), while in V^^ all the 
elements are equal to the corresponding elements of V except in the two central 
rows (n =^ ± O). In the row n = all the elements are zero except the central two, 
which are vi{m, = O) and V2{m = — O). In the row n = —0 the general element 



IS C^_o Cj^Q. 



§ 6. The Complete SoLUTioisr for Maxwellian Molecules. 

In the case of Maxwellian molecules, ^.e., molecules which are point centres of force 
varying inversely as the fifth power of the distance, the solution reduces to finite 
terms. This arises from the fact (q/. § 9 (C), 'Phil. Trans. ,^ A, vol. 216, p. 323) that 
for such a law of inter-action the functions 9^, ^, x ^^® independent of ?/, with the 
result that for all values of t we have {of, (4*25)) 



(o 01) /Cf — 1, >fcii — ^ iij ^ 12 -" ^ 12? fc 22 — ^ 

and consequently, 

s s 

(6*02) ll,^Gf^j^~*^2Jct — ^,Gt^i"*^2^ ™ (^1+^2)^ — 1- 




22? 







From (5*12)-(514) and (4*27), (4*29), with the aid of (6*02), we deduce that 
(6*03) aoo ^ a^s = ^ro = 2T -^^—^ \\2^oK'K-\2 (0). 

independently of the value of r or 5. From (5*21) and (3*13) we consequently have 



^'0 .. „ ^'0 .. _ ^'0 



OLq — , CL_Q 



(6*04) a^ = (r 7^ O), olo^gl-o = — == "" x x ^ -0 - . , --0 - . 

Again, from (6*03) and (5'16), (5*17) we conclude that h^ — — h, for all values 
of ±r, ±s from to co. Hence, by (5*36), or more simply from (5'2l) — since 
in Vo (&^^) all the elements of the central column are zero, instead of all except one, 
the central element, as in V — we have 

(6*05) ^'0 == 0. 
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Also in V (fe^^) and V^ (&^n) ^^ ^^^7 i-^ this case omit the central row and column, all 
the elements of which are zero {if r 7^ O) except the central element, which is the 
same in each (6) ; we suppose that this is done, so that V and V^ become determinants 
which possevss no single central row or column, just as with V (c^„). We will consider 
the two central columns of V (6^^) and V (c^^)i corresponding to m = ± 0. We have, 
by (512H518), (4-35H4-38), and (6*03), . 



(6-06) s>0, ^»^"= ^Xk743(^i.'+i~*o.-+i) 



(6-07 



{6-08 



(6-09 



(610 



6-11 



(6-12 



(6-13 



s < 0, 6_o. = - II^BoK',2 (0) {v,+ 5mi.i (3^/+mi^+|mimA/)}, 
s > 0, fe_o, = - ^lil^BoK', (0) mim/ (1 +-R/), 

S < 0, &o, = T^T^ BoK'12 ( ) /X1V2 ( 1 + i^^ia), 

S > 0, Co, = II ^CoK',2 (0) [3,,kn' + 2M2»'2 ( 1 Omi + 3^2^12°)], 

S > 0, C_o, = ||27rAC2.'2CoK'i2 (0) ( - IOmi + 3MA2"), 

00 
S < 0, Co, = — j 27r^iI/iCoK'i2 (0) ( - 10^2+ 3m2^I2°), 

5 < 0, C„o. = ||tCoK',2 (0) [3.A2''+2Mi«'i {10^,+ SfJi^kM 

15 



All these quantities, it will be noticed, are independent of s; thus we see that in 
V (b^^) and V (c^^) the elements of the two centre columns take only four distinct 
values, all the elements of these columns which are in the same quadrant having 

the same value. 

By using the method of differencing, by rows it is easy to prove, as a consequence 
of the property of V (b^^) and V (c^^) which we have just established, that 

(614) ^, = (r>2,r<-2), y, = (r > 1, r <--l), ^ 
while (0/ (6-05), (614), (314)), • 

(615) ro ^ ^^00-^^-0-0 , ^_^ ^ ^^^"^;-^; r 



b_Q_o—b_QQ Q __ ___g _ \bQQ—bi 



(6 16) /3o = -/3i ^ T — 7— — 7 7 ,-' P-o - -P-»i - 
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Consequently, in the case of Maxwellian molecules, F (c^) and H (c^) each reduce to 
their first terms, while G{c^) reduces to two terms only, the solution being quite 
finite. 

§ 7. The Genebal Soltttion when m, : m, is Very' Labge. 

When the mass rrij of the heavier mplecules is very large compared with that of 
the lighter molecules, so that m^jnii niay be supposed zero, an exact general solution 
may be obtained in simple terms. This was first proved by Lorentz,^ in connection 
with the theory of electrons. His method is much simpler than that of this paper, 
from which, however, his results may, with little difficulty, be deduced as a very 
special case. The deduction will be described in some detail, since the knowledge 
of the exact solution throws valuable light on the convergence of our successive 
approximations ; it also leads to an expression for tt which is of interest, as being in a 
form which, so far as my knowledge goes, is new. 

We suppose that the effect of collisions between the molecules m^ is negligible, so 
that P22 and P22 may be omitted from our calculations ; if the molecules are rigid 
elastic spheres of radii g-^ : cr^, this amounts to the neglect of crila-i {cf. § 9 (/) ) or, in 
general, to ^22/^11- 

It may readily be seen from (413), (414) that if m^ = the value of W{m, n), 
with whatever suffix, is zero except when ^ = 0, and that 

..^^. rB^,,,(m,r^) = (2^^)-+^ BV2 (m, ri) = 2-^y-x^ 

1 B^2i2i (m, n) = {2y'Y^\ B\,,^ (m, n) = 2^^-x'^f\ 

We also require B^ (m, n) to the first order in w^, as follows : — • 

(7'02) X-^''''" ^^' ^^ ^ |mn^22"^+ V^"^"^'^- V, B^2U2 (m, n) - -|mriM22^^V"^-'% 
\'B\,,2i (w, n) -= |mr2./x22^+V^^^+'^-^^cc^, B^^^i (m, n) = |m72,/>t,2^+VV- 

From these, by means of (4*07)~-(4M2) we deduce the following expressions for 
P (r, 5), p (r, s), retaining only the terms of highest order : — 

(7-03) PrAriSi) = |-7r.,.,M22^-^^(2ra + r + s + f)(r-h^ + i-),^,K\,(0), 
(7-04) p,, {r,s,) = ^i,rv,v,f,,2^^^ (^+f ).+! (^ + f )AK^i2 (0), 
(7-05) p,, {r,s,) = - |7r.,.,M22^+^ (r +f ),^, (^+f )AK^,, (O), 

^ LORENTZ, ^Archives N^erlandaises/ 10, p. 336, 1905 ; ^ Theory of Electrons,' p. 268. A more 
general form of Lorentz^s theory is given in Jeans^ ' Dynamical Theory of Gases/ 2nd ed., §§ 314, et seq. 
I am. indebted to Mr. Jeans for pointing out the interest of a detailed comparison of Lorentz's theory 
with this special ease of my own, and in consequence I have rewritten § 7 with greater fullness than at 
first. 
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(7-06) P2i{r^2) = i^yiy2f^i2'-*'{r+s+iX^.^,K^.K\,{0), 

(7-07) /isM) = fi7r.i.,A',2'-+'(r5 + r + s + f){r + s + |),^.^,K'i,(0), 

■ /i2 (riS:) = ^5'rn^2M,2'+Vs{r+|),^2 {s + i)A^iK\2 (O), r>0,s>0. 

(7-08) -l MO,S,) = f|7r.,.2M22' (s + fX^, (--^o + T^^i/) K'l, (0) = A (SiO,), 

(7 "09) /,, (n.9,) = ^57r.i.,Mi2^*Vs (r + i), (s + |-X+A+iK',;(0), r > 0, s > 0, 

(7-10) A (r,^^) = f|^.ii^./xi2^-^^ (^+^+5)^^^^^^^^r+.+i g/^^ (0)^ 

For the particular values of r and ^ there dealt with, the expressions (4'26)-(4*34) 
reduce to the above in the present special case, as may readily be verified. 

We now substitute from (7*03)--(7*10) in (512H5'20), retaining only the terms of 
the highest order, and neglecting quantities of order m/ : m{^ ; the result is as 
follows : — - 

(7*11) 5>0, ao, = ^V^iVoWiM2AoK'i2{0) = ^oo, ^<0, ^o, = ^TrXiXsi^o^iMaAo^^K^iaCO), 

(712) r > 0, 5 > 0, a„ = — 'N„pn (r^), 

(713) r < 0, 5 < 0, a„ = -^VXi^^'-o^aMiAo /^"^g w5h\ K.^'n (O), 

714) (r < 0, s > 0), a„ = — 2V^i^2>'o»»2MiAo^rK'i2 (O), 

715) (r> 0, s < 0), a„ = ^V?^,X2»'o"ii/^2Ao^,K'i2(0), 

716) (r > 0), b,= 0, b== ^aoo, ^'r = -^\~Wv,m,f,,Bo^-^;^K'M, (r< -O). 

717) (s > 0), &. =0, 6. = - ^^^ gXiVo^'iM2Bo^J^ K\, (O), (s < -O), 

718) (r<-0, «>0), b„ = 0, (r>0, s<-0), 

719) (r > 0, s > 0), ^" = X ^5 ^A2«'o"i2MiBo / ^^Wg^^) 






Since, by (714), (715), a„(r <0, «> O) is independent of s, and a„(r> 0, s < O) 
is independent of r, it is readily evident that in the alternate quadrants of the 
determinants V {§„„a^) and V (Ca„„) all the elements vanish. Hence each V becomes 
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equal to the product of the two remaining quadrants ; in the positive quadrant the 
determinants are the same (r> 0, s> O). Consequently in this case (5 '35) becomes 

(7 •20) a'o = -XiXaS (a^-a^^) = - ^^»»--_.^™™^_^ _o 

27ri/oW2AoKi2(0) D 

where D is the determinant whose general element is 

(7-21) §^^a.^ = S.^±^^+^k^^^, r>0,s>0, 

and Do is its principal minor. 

Similarly, it follows from (5*36) and (716)-(719) that ^^o can be likewise expressed 
as the quotient by an infinite determinant (covering only a quadrant of the infinite 
plane) of its second minor, as follows : — 

(7 -22) ^', ^ - ^^^^' ^'^ 



27n'oW2BoK'l2 (O) D' 

Here D^ is the determinant whose general element <i„ is given by 

(7-23) (r > 0, . > 0) d„ = ^._,,_, 1 {-i^±f+iW k,^.-~K-k,+ 1 

rsl(r + f),(s+f). 



rs " {■^' + f)r(s+f), """ 
(7-24) r > d^ = d,,^ 4_i - a,-l) = - SX doo = 1. 

Thus it is the same as the determinant D of (7'20) except that the r*^ row and 
column (f>l) are divided throughout by r. The determinant D'l is the second 
minor of D^, i.e,^ the minor of the next element to the centre in the first row or 
column. 

Again, from (5 '37), we may prove that 

where Di is the second minor of D just as J)\ is of D^ It is easy to see that^ 
(7*26) ^ = ^. 



Likewise, from (5*39), we have 



00 



(7-27) . S (.A + .'2/8_,) 



27.0 D"i 



1 27rvim3oK'i2 (0) D" * 

* Cf. the footnote on p. 154, indicating that Sr (//!«,.+ i'2a_,.) is, quite generally, a mere multiple of /S'o. 



MON ATOMIC GAS: DIFFUSION, VISCOSITY, AND THERMAL CONDUCTION. 149 

where D" is the determinant which is identical with D' except in its first two rows 
(for which cf. (5 "38)), while D"i is the minor of the second element in the first >ot<) 
of D". 

§ 8. The General Solution foe the Case op Similar Molecules. 

Another instructive and specially simple particular case is that of a gas composed 
of two sets of molecules whose mechanical properties— mass and mode of inter-action 
during encounter—are identical. In this case all the symbols which we have 
distinguished by the suffixes 1 or 2, to indicate reference to one or other molecular 
group, now have the same value for either suffix, with the exception of v^^ v^ or Xj, Xg, 
which denote the numbers or proportions of the two molecular groups. Moreover 
(c/. (4'0l)-{4*04), (4'07)~-(412), or, for a simpler and more general explanation, ' Phil. 
Trans.,' A, vol 216, § 7 (H), p. 309), we have 

1 1 1 

(8 '01) -3 pil (n^l) = —^ P22 (^^2) = -2 P ('^% ^)^ 

n ^2 ^0 

(8 -02) pu (n^i) +P12 {r2Si) =- P21 (^2^2) +P21 (^1^2) = X1X2P (r, s) pr^ (n^\) = P21 (n^2), 

1 1 1 

(8-03) -1 pn (rx5i) = -5 p'22 (^2^2) == —2p (n s), 

Vi V2 Vq 

(8-04) p\2 (riSi) ■\-p\2 (n^i) = P21 (^2^2) + Ai (ns2) = \^2p (n s), 

where p{r, s) and p {r, s) are defined by (8*01) and (8*03). 

By means of these relations we may reduce the expressions for a'o ^^^ ^^0 to a 
much simpler form. Our operations may be performed on the actual determinants 
which express the general solution, but they could, of course, be equally well described 
as transformations of the general equations (5*08)-(5*ll), 

First considering V (5^,a^,), we add to the m^^ column on the right {m > O) the 
corresponding column (— m) on the left of the centre ; the new element on the right 
is now given by the equations 

(8*05) (m > 0, n > 0) S^^ (a«„ + a„^„) = - Ao^^„N,^„ {pn (nhn,) +/)i2 {m,n,)~\-p^ (^2^1)}, 

KLn^mn {M' + X^X^) p (w, u) 

= Xi ^ Ao4nN^„/> {m, n), 



^0 



Vq 



(8*06) (m > 0, s < 0) S^^ (a^n + <^-mn) = -^2 ^ Ao4»N««p (m, n). 



m 

Y 2 
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We next add X^/Xi times the n^^ row (n > O) to the { — nY^ row; the elements of 
the negative rows {n ^ O) thus acquire the following values : — 

(8-07) (m>0, n<0) ^..n j«.n + ^-m«+ ^'(a,^, + a_,_,)} - aoo - ■-'^ AJ^^oPvA^iOi), 

(8*08) (m<0, n<0) 4Ja^n+ ~^ (^m-nj = -^ Ao(J^,N^,pi,(n5i), 

These transformations do not alter the value of V ((5^«a^^), nor of V^ {S^n<^^n) when 
applied to this determinant, which continues to be identical with the principal minor 
of V in its new form. Since in each of these determinants the elements of one entire 
quadrant (m > 0, n <C O) are all zero, they may be expressed as the product of the 
two simply infinite determinants formed by the quadrants (m > 0, n> O) and 
(m <; 0, n<^ O). The former quadrant, however, is the same for V and V^, so that 
(5 '3 5) may be written (after a little reduction^) in the form 

where V^^ is the principal minor of V^, the general elements a^^^ of which is given by 
the equations 

(8-10^ a' := r? ^ 27pi2 (^> n) 

^ ' ""' "^"-i^.2---^(m + fU,(n-ff)„,,K^,(0)' 

so that 



(8-11) a'oo - 1 aU = ^mo2-'« 2 „C,^, = a' 



m 







Here, since mi ^ m^, the equation (4*08) takes the special form 

2(n^i) 



(812) ^ 

V1V2 



f4-l,S + l 



^h^W e--''-^'xy E li^" (y) {B* (r + 1 , s) + B^- (r, s + 1 )} + 2t/V* (p) B' {r, s)] dx dy, 

whence it is clear that V^, and consequently also a^o> is independent of the ratio vi : ^2 
{i,e., vi and V2 occur only in the form 1^1 + 1^2 or vq). 

It may readily be verified that a corresponding calculation in the case of V (S^J)^^) 
leads to the result 

(813) /3', = 0. 

In the same way we may determine the particular forms assumed by (5*37), (5'39), 
and (5 '40) in this special case. 

^ This reduction chiefly consists in cancelling out factors common to all the elements of corresponding 
rows of Vo and Vq. 
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§ 9. FiBST AND Second Approximations in the General Case. 

(a) General Remarks, 

In the general case there are no such simplifications as have been described in §§ 7-9, 
and for practical numerical purposes we have to be content with successive approxi- 
mations to the complete solution. The approximate formulae become increasingly 
complicated with each successive step, however, so that it is fortunate that a first 
or second approximation generally gives a close approach to accuracy. In the theory 
of viscosity and conduction in a simple gas ('Phil. Trans.,' A, vol. 216, § 11, p. 334) 
it was found that a first approximation gave a result not more than 2 or 3 per cent, 
too small, while the error after a second approximation was negligible. In the 
present theory of a composite gas the error of the first approximation may be much 
larger (up to 13 pe^r cent.) in extreme cases where the masses, densities, ol* diameters 
of the two sets of molecules differ widely (c/. §7). Such a case is worked out to 
a fifth approximation in § 13 (e), and as this is one in which an exact solution 
is possible by another method, the results there obtained throw much light on the 
general character of the convergence of our analysis. It would seem that (as in the 
previous memoir just cited) the successive approximations form a monatomic sequence, 
the first and second members of which give a good indication of the accurate limiting 
value. If the difference between these two members is about 2 per cent., the 
additional correction due to all further approximations is about ^ per cent., while 
if the difference is so much as 8 per cent., the further correction is about 4 per cent., 
the additional correction being in an increasing ratio to the first difference as the 
value of the latter rises. 

In this paper we shall not go beyond a first approximation in the general case, 
except in regard to a'o? where we shall stop at the second approximation. 

(6) First Approximation to aV 

If in our set of equations (5*08) we neglect all save the central one, and consider 
only the central term of that, we get the equation 

(9'0l) aoo(ao-"a„o) ^ 1, or aoo^V^ -\K 

By (512) and (4*27) we have 
(9-02) a,, := Ao 2^5^^^ 

It is convenient to choose Aq (which is quite arbitrary as yet) so as to make 

(9*021) aoo-1, 

so that 

(9-022)' A,= ^7im, + m,) 

27rXiX2VQmim2isJi2\0) 
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and so that the equation for the first approximation to a\ becomes (with this value 
of Ao) 



(9*03) 



a 



AiA^. 



(c) Second Approximation to a\. 

We may conveniently obtain our second approximation by means of (5*35), taking 
the three rows and columns of V and V which contain aoo in the centre. From 
(5*12)-(5*14) it is easy to see that the elements have the following values (when 



(9*04) 4aoi = /^^(^j-l) = 4^10 =i>i, 4-i%~i = Ml (^1-1) = -~(5_ioa_ 



10 



■P2, 



(9*05) 



^ii«ii = if m/ + MaM 1 ■- 2^1 + ps) + -2%Mi/i2 






'12 1" "21" " '^'ii 



P 



Il> 



(9-06) a_z_ia_,_, = -^ jlW + Mi^(l-- 24 + 14)+ 2V1M2V + A -^ ^24 



"P32> 



(9*07) 



c5l_iai„i = -<f_na„„ =- MiM2{~V— -24+-p2~-2%42'^} =r -pi2. 



Consequently, the third-order determinant for V is given by 



(9*08) 



i>i2 Pi Pn 



P2 1 Pi 



J^22 P2 P12 



{Pu-PllP22'-2PlP2Pl2+PllP2''^P22Pl^} = ^Aj, 



and its principal minor by 
(9'081) 



P12 i^ii 
P22 P12 



{Pi2-PnP'i2)- 



It is easy to prove that 



(9-082) 



where 



1 1 






(9-09) d, = - hn' {Ji^(1-2m:/«.) + Mm,M2^,/-2 (H2^i-|4) Mi'}, 

Ms 

(9-10) d, = - V {i^^{l-2MiM2)+M-«iM2/fci/-2(i + 2A;i-P,);u/}, 

Ml 



(9-11) 2c^ 



— 4 
12 == 1"'5 



M1M2 



^11' V +1(1- 4miM2) ( 1 - 2^, + p,) + 8/x,M2 V (^5^ -- 2k, + P,). 
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Also 
(9-12) 



where 






251^1^2 



(913) 61 = 2/xi/ii2^ 



11 ) 



^2 = 2/X2M21^23^ 26j2 ^ 30(l— 4/Xi/i2) + 8/Xi/X2^,; 



Consequently, if we write 



(914) 



eo = 



2P1P2P12 -Pii P2' - P22P] 

J^12 '~'PllP22 



{k,-iy- 



diVi + 2(ii2^1*^2 + C?2»^2 



we have, as our second approximation to a^o by (5*35), 



{9-15) 



ao — "~X|X2 



1- 



We may note that 



(916) 



Ai = --i" {c?il^i'+2(^i2^il^2 + C^2^2')-(^l-l)M&i^l' + 26i2»^ll/2 + &2»^2')} 



V1V2 



1 



V1V2 



where 
(9161) 



(9162) 
(9163) 



o 

a, = —kn' {|m2'"-Mi'(V— p2) + ^VlM2^12'}, 
M2 

0^2 = ^ ^22' {fMl'-M2M^l'~-|^2)+^%/tXi/X2^12'}, 

Ml 

^12 ^^ ^12 1^1 1/ ^12» 



The determinant Vo(^^^a^^) of (5 '37) is, to the same order of approximation as for 

V(^^,aJin(9-08), 



(9164) ~ 



Pu Pi Pn 



V2 Vi 



P22 P2 Pu 



^1 {P2P12-P1P22) + V2 {piPi2-'P2Pn) = -2^5^ (^1-- 1) ^1 



where 



(917) €1 = 30 (j'iM2 + '^2Ml) (mi-M2)-8/XiM2^^2 {Vl-i'2) -^^ {^lM'12f^\l-^2^2lk%2)' 

Hence we have, as a first approximation, 

7 1 

(9-18) ir{v,a, + ,,a_;)=^e,. 

1 ^i 



154 



DR, S. CHAPMAN ON THE KINETIC THEORY OF A COMPOSITE 



{d) First Approximation to ^^q and 2{i/i^^— i/^jS^^). 

In this case, owing to the fact that the central element in the numerator 

determinants of (5*36) and (5*39) is zero, we hare to deal with the third order 

determinant in order to obtain a ^st approximation. From (515)-(518) it is easy 

to see that V {^^„6,„„) may be obtained from V {S^n<^„,n), each being limited to the 

three central rows and columns, if we replace c^oo by x-%<^oo5 and divide the top and 
bottom rows respectively by Xi and — Xj, Consequently the third order determinant 
obtained from ^ {S^nKm) is equal to -- iV— (t-w) ^i- The value of V^o(4«&in») is 
similarly {cf. § 5 (/) for the definition of this determinant), 



(919) 



B, 



.-O-aXw. 



Pl2 1 £1.1 
Xi Xi 



P2 



23 



X. 






Pi 

^2 



1 



XiX2l'o X-A-oii/ 



P12 "1 Pn i 



7 1 / T> 



2 







P2 p, 



J^"^ £/& jS ,.£/ 1 -d 



25X1X2% \A0Ry 



€3 



by (917). Hence, by (5*36), we have 



/3'^ 



JBa V 



0^1 



-D 



'0 "0 1 



(9-20) j3\ = - ^~ -^^i-^ XiX,e3 

by (918). 

Again, the determinant V(, (<^„„6'„„) defined in § 5 {g) has the form 



(9-21) 



JJq 



, juLaJLv/ 



L Pv. 



X, 



■P2 



1 Pj2 



— -^0 






iXXaXw' 



^2 

Pm 

P2 



'P2 



1 21^P2 







.I/Pl^pi 

\Xi X2. 



\ 



1 



1 



2\x .1 / Pl 



^1 



^1 Ag 



1 







Pl 



V-i 






Xi X2 



vo 



'0 



^1^2 \-^^0 



XJu-aJlv 



Pl^^222 ^ 
\ ^2 

2 

' iP22 -Pu') + 2p,v, {P12-P1P2) + vi (pn -Pi')] 



yl 



vo 



B 



'0 



-%1 ^22 



C^2^^ 



This last equation, here proved true only by comparing the first approximations to ^'q and 
^•r{n^r+n^~-r), may easily be shown to be strictly accurate, by comparing the general expressions for 
these quantities. 
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where ai and 0^2 ^^^ given by (9*161), (9*162) and a\2 by 

(9-22) a\, - M1M2 jf->t/+ p2-^%^i/ + ^B (^+M 

Consequently, by (5*39), to this order of approximation we have 
(9-23) i(.,^. + .,/3_,) = ||^{^oa,.,»+2aV,.,+ f^^a,./' 



2B0V0V2 1^11 k 





22 



(e) JT^e i^^rs^ Approximation to 2(viy^ + v2y-r)- 
The two central terms of the two central equations of type (5 '11) are 

(9-24) y^\l+^^f,^^k,U^o'^—k\}--y__,{l-^ak,')=J- ^^TrFT-fTTy 

(9-25) -yo(l-^A/) + y-ojl+^aMi2-fci/+2%^-^^24 = ^ ^WW^' 

I »^1 M1M2 J 07n^il/2/iiM2H-^ 12 \^/ 

which lead to the result (equivalent to that obtained as a first approximation 
from 5*40), 

^y ^bj Oo 2. l.,y.4-.2y-.; - 2^g^^2 (0) ^11%.1^ + 2"&12V,2.:.2 + ^^22^2^' 

where 

(9-27), {9-28) Ci = 1 -f^%MiA/, c, = 1 +-^^^2^/, 

(9-29) c,,= l-^^h,' + ~is-^{k,' + hj), 

(9-30) c\,= H-^-LMl. 

/X2/X2 '^]2 

A second approximation would replace the right-hand factor in (9*26), by the 
quotient of one quartic homogeneous polynomial in v^, V2 by another, and similarly each 
further approximation would increase the degree of both numerator and denominator 
by 2. 

(/) Values of K^i2(0) for Particular Molecular Models. 

In connection with the above approximate formulae, and their subsequent 
applications, it is convenient here to write down the values taken by K^i2(0) and the 
other constants k, when certain particular characteristics are attributed to the 
molecules. 
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When the molecules are regarded as rigid elastic spheres of radii a-^ and 0-2, it may 
readily be shown (cf, * Phil. Trans./ A, vol. 216, § 9, p. 320) that 



(9-31) K\,{0)^i{<T,-^a-,y 



\h7rmim2/ 



When the molecules are point centres of force proportional to the inverse n*^ 
power of the distance, the force at unit distance being Kj^mi^n^, it may be shown"^ 
that 

(9-33) K'.,(0) = 4,1. W(|!5^j*""'[K„(™,+™.)]A r(3-^- 

where 

(9-34) ^ Ix(n)~47r mn^e.ada 

Jo 

and 



(9-35) ^^0^2\'—~rc--^^ 

1--^^---"— r(~) r 

n— IW J 







VI 



>;o being the root of the following equation in tj : 
(9-36) l-^^-~^f^Y'' = 0. 



Further, 



n—l \a, 



T(t + 3 ~] rlt + i ^ ^' 



(9-37) k = ^ , ^-^^^ >^,/ = 5 - - ^- Z^-^V^^ I. (^) 



where 



(,+3Xr(3-^) (*+IX.xr(3-^^i^(«) 



QO 



(9-38) l^in) = TT sin' e. a da, 

Also 



*'0 



2 2 

K\n~l /TT \^-i 



W ^"/ ^11 — (iF~~/ '^12 ^^22 — ("^"^i ^12 J 



where Kn, K22 are the force constants appropriate to a pair of molecules of the first 
or second kind respectively ; we here assume that the law of variation of force with 
distance is the same whether the molecules are like or unlike — if this be not so, n must 

* Cf. ' Phil. Trans./ A, vol. 216, § 9, p. 320, or Jeans' ' Dynamical Theory of Gases ' (2nd edit.), §§ 305 
et seq, I have adopted the notation of the above equations in order to facilitate comparison with the 
corresponding equations in Jeans' treatise. 
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also have suffixes 11, 12, or 22, and the formulae (9;37)~(9*39) will need to be modified 
in a manner sufficiently obvious. 

When n — 5, i^e.^ when the molecules are Maxwellian, it is clear that 

(9-40) K',,{0)=|-I,(5)[K,,K+mJ^ 

(9 41) Aj ™ 1 k^2 " 5 -^ , V . 



§ 10. The Equation of Diffusion. 

(a) Definitions of the Various Coefficients of Diffusion, 

In our analysis of molecular motions (§1) the rate of inter-diffusion was expressed 
in terms of u^q, which is defined by the equations 

(lO'Ol) (j^l + V2)tt^o = ^li'^l—'^o) = •"^2(^2-'^o) = -"^^^^(^1~~'?^2)- 

Vi + V2 

Further, by (315) and (2*06), we have, as the equation of diffusion, 

(10-02) „.. = jA..'.(±|^.-X'.-^^)+iBx||. 

Hence it appears that the agents effective in causing diffusion are (taking the 
terms of this equation in order) a concentration-gradient or variation in the relative 
proportions of the two component gases, external forces of unequal amounts per unit 
mass of the two gases, and variations in the total pressure or the temperature of the 
composite gas. 

When the pressure and temperature are uniform, and the external forces are such 
as to make X^o> Y^o? ^^o zero, diffusion can take place only if there is a variation in the 
relative concentration of the component gases. In this case we may compare (10*02) 
which now has the special form 

(10-03) < - iA,a^o ^ -^ -= TT- ^0^ oT HT = "- IZ~ ^0^ 51 ^ 



2h dx 9.0 ' '2hdx 9^0 2hdx' 



with the ordinary equation of diffusion 

■J 
(10'04) vqU^q = viiui-u^) = -Di2-^ v^iu^-u^ = -D 



I vAu,-u,)=-dJ^, 



where Djg is the coefficient of diffusion. By comparison we have 

(10-05) D,, = --iAoRTaV 

z 2 
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When the pressure, temperature, and composition of the gas are uniform, a steady- 
motion of diffusion may be produced by equal and opposite forces ± {voX!q, j^qY^o? ^o^^o) 
per unit volume on the two components. This will be termed forced diffusion, and 
we define the coefiicient of forced diffusion D^g by the equation " 

(10-06) ^^-DVoX'o, 

so that, by comparison with (10*02), 

(10*07) I)\, == ~~i i A,a^o - - Di2. 

From (10*02) it is clear that diffusion will occur also when the relative proportion 
of the two gaseous components is uniform, and in the absence of external forces such 
as might produce diffusion, provided that the pressure or temperature varies. If we 
define coefficients of thermal diffusion Drj. and of pressure diffusion D^ by the 
equation 

(10-08) «'^==D,i-|^»-D4|^, 

pQ OX i ox 

we have, by comparison with ( 10*02) 

(10*09) D, - ^|Aoa;£^^- - -iAoETa^o ^ - ^Di2. 

romo mo mo 

(1010) D, - -iBo/3'oT. 

(b) The Equation of Diffusion. 

If we now substitute the various coefl&cients of diffusion in the equation of diffusion 
(10*02), this becomes 

(1011) «'^=._D,|^» + DUX'o+D,i|^»-D,iF- 

OX ^ Pq dx 1 OX 

In later sections, when we consider in detail the values of the four coefficients of 
diffusion, we shall see that they are all positive (the molecules 1 being the heavier — 
cf. §1)."^ Hence from (10*1 1) we deduce that the direction of diffusion of the heavier 
component of a gas is (a) opposite to the direction of increasing concentration (6) in 
the direction of the diffusion-component of the external force (c) in the direction of 
increasing mean pressure and (d) opposite to the direction of increasing temperature. "^ 
This is fairly evident as regards (a) and (6). In case (c), the sign of D^ is the same 
as that of m^o (^- (10*09), Di^ being always positive), which is a multiple of mi— m^; 
the physical reason is also not difficult to grasp— under the influence of a difference 
of pressure both components will tend to flow in the direction of diminishing 
pressure, so as to render the pressure uniform. The lighter molecules will travel 

"^ This statement is modified, with regard to D^, in Note E, p. 197. 
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faster than the heavier, however, so that relative to the mean motion the heavier gas 
will be diffusing in the direction of increasing pressure. Such a process must go on, 
to some extent, during the passage of sound-waves in air or any mixed gas, and will 
influence the propagation of the disturbance to a degree which is probably comparable 
with that of the effects due to viscosity and thermal conduction. The effects of these 
latter have been examined by Stokes, Kirchhoff, and Eayleigh.* 

It is rather more difficult to perceive, either analytically or physically, that D^ 
must be positive if the molecules 1 are the heavier (see Note E, p. 197). A discussion 
of this coefficient will be found in § 14. 

If in (10*11 ) we assign to u\ the value zero, we get the equation of condition for a 
state without diffusion. Some special cases of this equation will be considered in § 16, 
after the relative magnitudes of the various coefficients of diffusion have been 
determined. 



§11. The Equation of Viscosity. 

The various pressure components p^^^ p^j^ et cetera, for the composite gas are given 
by the following and similar equations : — 



(iroi) 



Pxx = ^"^1 U / + v^m^ U 2^ Pxy = »^i^i U 1 Vi + v^m^ U 2V2. 



A reference to (5'24) or to our expressions for /(U, V, W) in § 3 hence enables us to 
deduce that 



(11-02) 



Pxx-P 







15 



1 ^ 



2h 







(iro3) 



p 



xy 



¥ 



1 ^ 

7 5 ^ Co^:cy S (l^iVr + J^sV-r)? 



We may compare these with the equations of pressure of a gas whose coefficient of 
viscosity is ^12(0/. (2*11) and (2*14) for the values of c^^ and c^^) : — 



(11'04) 



(iro5) 



Pxx'^Po -~ ~^^12 



dx ^\dx dy 



32 



— Sf^U^xxf 



P 



'dVa . du 



xy 



'^12 



_|_ _ ] _ '3,tCi2Cy.y, 



dx dy 



It thus becomes evident that the composite gas behaves like a viscous fluid whose 
coefficient of viscosity is given by 



* Stokes, 'Cambridge Transactions,* 8, p. 287, 1845; Kirchhoff, *Pogg. Ann.,' 134, p. 177, 1868; 
Eayleigh, ' Phil. Trans.,' 175, p. 1883, * Theory of Sound,' II., ch, XIX. 
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§ 12. The Equation of Energy. 
(a) The Deduction of the Equation. 

In order to obtain the equation of energy we return to the consideration of the 
equation of transfer (2*01). Having determined the form of the velocity distribution 
function/ (U, V, W) correct to the first order of small quantities, we can now make 
the equations of transfer accurate so far as the second order, and this is necessary in 
connection with the equation of energy. Also we need no longer suppose that the 
mean motion of the gas {uq, Vq, Wq) at the point considered is zero. 

We add to (2*01) the corresponding equation for the molecules of the second kind, 
with the result 



p^ r" /^ 

(12*01) ■^^Mi + y2Q2)+^ g^{^i(^)iQi + .,(^),Q2} 



Ml \S{u)J ma \3(ti) 



2/ J - 



= A(Q, + Q,). 



In this equation we substitute in turn Q = 1, Q = m (?/), and Q = ^m {{uf -{- {vY -\- {wf}, 
thus obtaining the equations of density, momentum, and energy. Since all three of 
these quantities are conserved unchanged during the molecular encounters, we have 
in each case A (Qi + Q2) = 0. The first equation takes the form 

/12-02) ^ + ^^^ + ^^ + ^-^^^ = 

where we have made use of the notation of § 1. We will denote the operator 



a 

dt ' ^^"^ dx ' ^""dy ' ^"""dz 



— + t^o;r^ ■^Vq — -\'Wq 



D 

by Y^ y commonly known as the '' mobile operator " ; if q is any function of 

^) 2/ ^? ^j ^ denotes the rate of change, with respect to time, if the value of q at 
Dt 

a point which moves with the mean velocity of the gas. Then the equation of 
continuity (12*02) may also be written in the form 

(12-03) lDilo^/^+|^^.^\^0. 

•^ ^ vq Dt \ ox dy dz I 

The second equation, the equation of momentum, may be written as 

>\ '^ ^ ■ p\ 

(12*04) - {piUi-\- p.u^) + -^ [pi W-Vp^iu),'} + g- {pi {u)i{v,)-\-p, {u), {v,)} 



dz 



MONATOMIC GAS : DIFFUSION, VISCOSITY, AND THERMAL CONDUCTION. 1 6l 
or (c/. § 1 and (ll'Ol)) as 

(12-05) l/flXo = ^ (/)oWo + p'oO + g^ {/5oMo^+ 2p>oW'o+Pxx} 

dy 

\ 3a; 3?/ 3^; / \ dx - dy dz 

In deducing the third equation, the equation of energy, it is convenient, first of all, 
to write out the following equations giving the appropriate special values of the 
various quantities which occur in (l2-0l) :- 



(12-06) viQi + i/A = h^i{e^ + ^—^u^o-^G^) +ii^2^2K^-2^ St^o^+C/ 



(12*07) vi {u)i Qi+ ^2 M2 Q2 = "^Q HpqOo^+poS.UqU^q-^^'RvqT) +i/o'^Vi 



2 





+ {^oPxx + %Pxy + '^oPxz) + (i/>lUiCi^ + h^^f^i)^ 



(12-08) J^Xx(^-V) + -^X,(^) = v,u,lL,^v,u,X, = v^oXo+Vo<(X,^X,). 
mi \3 \u)i/ mg \3 (t^ja/ 

On substitution of these values in (12*01), this becomes 
(12-09) V, ^ {hy^,G,'-v{m,-m,) 27>6o<+fRT} +^2 ^ (pVoVo) 

+ 2 — (ttoi>x:. + VoP^y + Woi>x.) + 2 ^ (i^i U iCi* + i^a U 3O/ ) 

2t6'o (Xi— X2) — vo-St^oXo = 0. 



Vq 



This equation can be simplified by eliminating vqZuqXqI thus from (12*05) we 
deduce that 

(12-10) voZu,Xo = vo^^{hhCo')+^o^^ojy^{^i-^n2)u'o + PoOo'^ 



\ 9a? 3^ d2J / ox 



2 
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by means of which we may transform (12'09) into 

(1211) fR.o §^ + vo (mi-m,) S«'o ^ -hW (^ + ^ + ^-") -"o^^ (X,-X,) 

In this equation we substitute the values of p^^., p^^, p^^ from (11*04), (11*05), et 
cetera. Also we divide throughout by J, the mechanical equivalent of heat, and 

replace f -^ — by C^, the specific heat of the gas at constant volume."^ After a little 
reduction we obtain the equation of energy in the following form :- — 

(&) The Interpretation of the Equation of Energy. 

We will now consider the significance of the various terms in the equation of 
energy (12*12). The terms on the right depend on difiusion, viscosity, and thermal 
conduction (as will be made clearer later). If we neglect the small changes of energy 
produced by these means, as a first approximation (12*12) may be written 

^^^^^^ HDt^Adx^ dy^ aj~^' 

where we have omitted the right-hand side of (12*12), and divided the left hand by 

* The specific heat of a simple gas at constant volume is f y— , if m is the molecular mass (c/. Jeans' 
' Dynamical Theory of Gases,' 2nd edit., § 261, (512)). Hence for a composite gas we have 

and consequently 
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By the equation of continuity (12'03), the last equation is equivalent to 

(12-14) ^S^-f-S^"=0, 

or to 

(12-15) ^^{T,,-'U) ^ ^(po^o-''') = 0. 

Since we are neglecting diffusion, mo (or Aimi + Xsmg) is invariable, and (12' 15) 
expresses the adiabatic law of expansion of a monatomic gas as a given element is 
followed throughout its motion, viz.. 



(1216) i>o = W^ T^Jcfpo 



% 



Thus (12*12) may be regarded as giving the correction to this law owing to 
diffusion, viscosity, and conduction. The left-hand side represents the net rate of 
increase of energy of molecular agitation, being the increase corresponding to the 
rising temperature after allowing for the change of heat energy produced by adiabatic 
expansion or compression. The right-hand side indicates that this is due to the 
following causes :— 

(a) We are considering the fluid contained in a volume element which moves with 
the mean motion of the fluid ; but in addition to this motion there is one of inter- 
diffusion within the element, there being a stream of ± {vqU^q, vqV^q^ pqW^o) molecules 
per unit time per unit volume in opposite directions. The cc-components of the 

effective forces on each of these molecules are Xj— mi -^-^ and X^— ms -jt^ respectively, 

the second term depending on -y~ being present, since we are throughout concerned 

Ut 

with the energy relative to the mean motion {uq^ Vq, Wq), The work done on the 

fluid by these forces is expressed, in thermal units, by the first right-hand term of 

(1212). 

(/3) Owing to the motion of diffusion there will be an increase in the number of 

molecules mj and an equal decrease in the number of molecules mg, at the rate of 

-TT-^ -4- ^r-^ + -TT"^ ) per unit time per unit volume. The mean energy of peculiar 

dx dy dz / 

motion is the same for molecules m^ and 7n2, so that this change affects only the 
energy of mean motion, which is -^rriiCQ for a molecule m^, and ^tti^Cq for a molecule 
m2. The increase of energy due to this cause is represented by the second right-hand 
term of (1212). 

(y) The third term gives the heat generated by the viscous forces acting in the 
gas, being, in fact, twice the '' dissipation function '^ of the viscous motion."^ 

* Cf, Lamb, * Hydrodynamics/ p. 518. 
VOL. OCX VII. — A. 2 A 
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(S) The fourth term represents the increase of heat which ordinary physics regards 
as due to conduction, together with a term due to diffusion which has not, so far 
as I know, been hitherto recognised in either ordinary physics or the kinetic theory. 
These two terms are discussed in detail in the next section of this chapter. 

(c) Conduction of Heat and the Thermal Flux of Diffusion. 
If there is no mass motion the equation of energy takes the form 

(12-17) poC„^ = i2.o<{X,-X,)-i 2 ^{hrU^G^^+h^um, 

and in the last term we will now substitute from (5*25), with the result 

9T 1 

(1218) PqG^ ~- =: j ZvqU^q (Xi — Xg) 



'^rn 00 



It is convenient to eliminate ^'q by means of (3*15), and on so doing (12 '18) 
becomes transformed into 

8T 1 
(12'19) poC^"^ = T Si/o^^o(Xi— Xs) 



+ 2 ART 

dx J 






'^^o + 9 Bo — \ 2 (i/a/3^ + ^2/3_^) 
oLq ox ll 

— i ^ 2 r {vxa^ + i/a^-r) 

If we suppose that no diffusion is taking place, so that {u^q^ v^q, w^q) are all zero, 
and compare (1219) with Fourier's equation of conduction of heat in a gas at 
rest, i,e,^ with 

we obtain the following equation for the coefficient of thermal conduction S- :— 

(12-21) ^=^iB,~\^{p,/3, + p,^^,)-^^r{via,-\-p,a„,)\. 

J L 1 oLq I J 

From (1219) we perceive also that the motion of inter- diffusion is accompanied 
by a flow of heat, which is proportional to the velocity of diffusion, and also depends 
on the temperature and the molecular densities, but is independent of the other 
characteristics of the gas. We shall term this process the thermal flux of diffusion ; 
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in the absence of thermal conduction and of external forces we shall have as the 
equation of energy 

(12-22) poa— = s|-(©o, 

where li is defined by the equation 



00 



RT 1 

(12*23) li = "~T~~T Sr{j/iGt^+j/2a„^X 

and is termed the specific energy of diffusion. 

We now proceed to consider in detail the various coefficients of diffusion, conduction, 
viscosity, and specific energy of diffusion for which we have obtained general formulae. 

§13. The Coefficient oe Diffusion Djs. 

{a) The General Formula. 

The general formula for the coefficient of diffiision Dj^ is obtainable in terms of the 
molecular data by substitution in (10'05) of the value of a^^ given by (5*35), thus 

(13-01) D,, = iAoETX,X,^|||^ - 

So far as § 9 (6) Ao had remained quite arbitrary, only the product of Ao into the 
coefficients a being definite. We there defined A^ as having such a value as to make 
aoo equal to unity. We now substitute that value {cf. (9*022)), viz., 

(13.02) A = 27(mi+m,) 

in (13*01), with the result that 

27n/o^im2K'i2(0) V(^^„a^J ' 
where v is the prmcipal minor 01 V. 

(b) The Case of Maxwellian Molecules. 

In the case of Maxwellian molecules it is readily evident from (6*03) that all the 
elements of the first row and column of V (4m^m»)5 save the initial element a^^^ are 
zero. Hence in this case V is equal to a^^V^ i.e., 

(13*04) T i^n.n(^rnn) J:^ ^ ^ 

V {O^j^ajnnj ^00 

2 A 2 
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in view of the convention in § 13 (a) as to the value of Gqq, Hence, by (13*03), we 
have 

3 (m|4-m2)RT 



(13-05) Di2 == 



27rvom,m^K\2{0) 



as the exact value of D^g when the molecules are Maxwellian. The same result 
follows readily also from (10*05), (6*03), and (6*04). It is the same, except as regards 
the notation, as the formula deduced by Maxwell in his second great memoir 
on the dynamical theory of gases. ^ 

(c) A First Approximation to the Coefficient of Diffusion, 

Only in the case just considered does our general formula for D^^ reduce to a 
simple finite form : usually we must, for practical purposes, be content to make 
approximations to the exact result. As in § 9, this may be effected by taking 
successive finite convergents to our infinite determinants, which is equivalent to 
neglecting all terms in the expansion of /(U, V, W) after the first one, two or more 
at the beginning. Thus for a first approximation, taking only the central element of 
"^ {^mn<^mn)j it Is cloar that VyV is equal to l/^oo or unity simply {cf, § 13a). Hence 
we have 

(13-06) D,, = i^^^;^^^ (1st approximation), 

a result which also follows from (10*05) and (9*03). This, it will be noticed, is the 
same as (13*05), showing that what is in general only a first approximation to D^g is 
in the case of Maxwellian molecules a strictly accurate result. 

The formula (13*06) is not new ; it was first given by Langevin,! and subsequently 
by myself I independently. In all these cases, and also in Maxwell's investigation 
the method used was an approximate one which involves the assumption that 
the peculiar velocities of the molecules of the two constituent gases are distributed 
about the separate mean velocities Cj, c^ according to Maxwell's law for the steady 
state of a gas. The method of the present paper is based on an actual determination 
of the law of distribution. The assumed law just mentioned, if expressed in the 
manner adopted in § 3, would involve no a-coefficient beyond a^ (neglecting squares 
or higher powers of the velocity of diffusion c^^). This coincides with the true law 
only in the case of Maxwellian molecules (§ 6), so that only in that case is (13*06) 
exact : Maxwell himself did not prove this rigorously, though he obtained the 

, * Maxwell, ' Phil Trans.,^ vol. 157 (1866) ; or ' Scientific Papers/ ii., p. 27. 

t Langevin, *Ann. de CMmie et de Physique' (8), v., 245, 1905: cf. also Enskog, * Phys. Zeit.,' xii. 
533, 1911. 

t ' Phil. Trans./ A, vol 211, p. 499 (35), 1911, 
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accurate expression. In the general case (13'06) is an approximation only, and the 
method used by Maxwell and the other writers just named does not afford a ready 
means of estimating the error involved. The present investigation makes it clear 
that in extreme cases this is far from small. 

(d) A Second Approximation to the Coefficient of Diffusion, 

Passing now to a second approximation, by means of (10*05) and (915) we obtain 
the result 

(13*07) Di2 = :; {m^-i-Trh) ^gnd approximation). 

^ ^ 1— 6o ^TTi/omirngKislO) 

The correction to (13*06) consists of a factor (l--6o)'"\ eo being given by (9*14) and 
(9*09)-(9*13) ; since for Maxwellian molecules A^i = 1, it follows from (9*14) that 
eo = in their case. In order to estimate the importance of cq in other cases we must 
consider some special typical molecular models, and determine the numerical values of 
eo for such gases. The most important models are those for which K\2(0), k^, etc., are 
given in § 9 (y). 

In paragraph (g) of this section we shall consider in detail the range in the values 
of eo calculated from these numbers, for various ratios of the molecular masses and 
diameters. Since cq only affords a second approximation to D12, however, and not the 
exact value, it is convenient first of all to examine certain particular cases of our 
formulae (c/! §§7, 8) which throw some light on the accuracy of a second 
approximation in general. 

(e) The Coefficient of Diffusion when mi : m2 and ^1/0-2 are Very Large, 

In the special case when the mass and size of a molecule of the first kind are so 
great, compared with that of a molecule of the second kind, that ma/mi is negHgible 
(c/. § 7), we obtain from (10*05), (7*20) and (13*03) the result 

/I ^-08 q\ n - ^^"^^ Do __ Do 3 

(13 08, 9) .U,2 - 2x.om2K\2(0) D ^ D IQv.ia.-Va^f {h^m^f^' 

where D is a determinant (occupying only a quadrant of the infinite plane) whose 
general element is given by (7*21), while Do is its principal minor. As usual, for 
MaxweUian molecules Do = D. We proceed to make successive numerical approxi- 
mations to Do/D in the case of molecules of other types. 

For molecules which are rigid elastic spheres we have seen (9*32) that 

/ X 7 (^+2), 

{13-10) *'=(Tfi-:' 
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so that by (7'2l) the general element of D is 



(1311) 






'r+s 



'rs^rs 



r>0, s>0. 



The following are the numerical values of the elements as far as the fourth row and 
column, the exact values being given in the first expression, and the decimal values 
to three significant figures in the second : — 



(1312) D = 



1 


1 
5 


1 

5.7 


3 
5.7.9 


— 


I'OOO 


0-200 


1 


13 
5.5 


23 
5.5.7 


33 




0-200 


0-520 


5 


5.5.7.9 


1 
5.7 


23 
5.5.7 


433 
5.5.7.7 


1077 
5.5.7.7.9 




-0-029 


0-131 


3 


33 


1077 
5.5.7.7.9 


26613 




0-010 


-0-021 


5.7.9 


5.5.7.9 


\J » \j » t » i * «/• <7 



0-200 -0-029 0-010 



0-021 



0-353 0-098 



0-098 0-268 



By neglecting the rows and columns of D after the first, second, third and fourth 
respectively we obtain the following four successive approximations to Do/D : — 





Table I. 








Approximation. 


D. 


Do/D. 


1st 


1-000 




1 • 0000 


833 


2nd 


13 

25 -^-^^^ 


13 
12 


= 1-0833 


235 


3rd 


4608 
25.25.49 - ^ 1^^ 


5100 
4608 


- 1-1068 


97 


4th 


106168100 

567492 — ^ ^^^ 


11854084 
1061681 


- 1-1165 





The successive approximations to Do/D evidently converge with some rapidity 
to the value 1*12 or, more nearly, 1*13; the correction introduced by the second 
approximation covers two-thirds of the error of the first. 

It is of importance and interest to notice that in this case an exact solution is 
obtainable by another method, less general than that of this memoir, but more 
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effective in the particular case before us ; this solution is due to Lorentz, and may 
be expressed in our notation by^ 

(13-13) Di2 = 



This is identical with our own solution (13*09) if 

(1314) §p = 32 ^ 11317,t 

^ ^ D 97r 

and our approximations to Do/D show how this numerical value is approached by the 
successive convergents to our infinite-determinantal solution. | Incidentally, we may 
notice that (1314) affords an expression for ir of a kind which, so far as I am aware, 
is quite new. 

If the molecules are n*^ power centres of force {cf. §9 (/)), our general formula 
(13 '09) becomes (when Wh2lwhi is negligible) 

(1315) D,2 = ^° ^'^ 



^ 8 i/oli {n) {hm^m^^^y-^ {h-irm.f' T ( 3 - ^ 



n—1. 



and the general term of the determinant D is given by 



r r+s+3 



(13-16) .5„a„ = 8„ 1 J" ^ 



(r+f)r(s+f).r(3- 2 



\ n—1, 

With this we may compare the exact formula obtained by Jeans {loc, cit), using 
the method of Lobentz, 



fr(2+ 

(1317) Br, = 



n—1 



k 



* Gf. LoRENTZ, ' Archives N^erlandaises/ 10, p. 336, 1905 ; /Theory of Electrons/ p. 268. Also, for a 
more general theory along the same lines, cf, Jeans' * Dynamical Theory of Gases ' (2nd ed.), § 333 (654) 
or § 450 (890). 

PiDDUCK, 'Proc. Lond. Math. Soc' (2), 15, p. 112, 1915, has also deduced Lorentz's result as a special 
case of the general method of solution by integral equations. 

t [I have now obtained a rigorous proof of this equation and of (13*18). — February 22, 1917.] 

I Cf. the footnote to page 171. 
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These are equivalent, so that 



J) 



^t(2^-^)tIs- 2 



'TT 



n—l, 



n— 1. 



16/l+^^_l_f2 2 



9 



TT 



n— 1/ n—1 



n—1 



i^^A^rl-^vd ^ 



,n— 1 



16/i + _2_^_2_^2-. 2 



W— 1/ 72-—! \ n—1. 



1- 



n— 1, 

2 



n—1/ • 27r 
sm 

n—1 



n—1 



by a well-known formula in the theory of gamma functions. When n = 5 {the case 

of Maxwellian molecules), the last equation gives Do/D = 1, as it should do. When 

32 
n = (x>, corresponding to the case of rigid elastic spheres, Dq/D = — , as before. 

It is of interest to consider one or two intermediate values of n in order to see with 
what rapidity our determinantal expression for DJD converges to the value given 
by (1318) ; we shall not go beyond a thiyd approximation. 

The determinant D, as far as the third row and column, is as follows, where ^m 

has been written in place of 



n—1 



(1319) D = 



1—m 



1— m 



13 — 4m+m^ 
5.5 



1— m^ 
5.7 

23-- 27m+ 5m^— m^ 



1-m^ 23-27m + 5m^-m^ 433-216m + 70m^-8m*+m- 



5.7 



5.5.7 



5.5.7.7 



In the following table are given the first three approximations to Dq/D for a few 
typical values of n, together with the exact values calculated from (13'18) : — 









Table II. 








n. 


m. 


Do/D. 




Approxi r r \ ations. 
















Exact value. 














1st. 


2nd. 




3rd. 




5 


1 


1-000 


1-000 




1-000 


1-000 


9 


i 


1-000 


1-023 




1-027 


1-031 


13 


1 


1-000 


1-039 




1-048 


1-056 


17 


i 


1-000 


1-049 




1-060 


1-072 


00 





1-000 


1-083 




1-107 


1-132 
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This table well illustrates the general remarks on the convergence"^ of our formulae 
which were made in § 9 (a). 

(/) The CoefficieMt of Self- Diffusion of a Gas. 

The next simplest case of a general character which we shall consider is that 
in which the two sets of molecules are alike in mass and dynamical properties ; 
a special case of this is that of a simple gas, the two sets of molecules being identical— 
the corresponding coefficient of diffusion D^g will be written Dn and termed the 
coefficient of self-diffusion. From (10'05), and (8'9), the latter equation being the 
one appropriate to the case before us, we find that 



(13'20) 



D 



3RT V^ 



11 



7n/omK^i2(0) Vq 



where we may now omit the distinguishing suffixes 1, 2; V^^ is the principal minor 
of Vq, whose general element is given by (8*10). We may note that the factor V^o/V^ 
is independent of i^ and 1/3 (§ 8). 

As far as the first two rows and columns, V^^ has the following form :— 



(13*21) 






1 



i(^l-l) 



2 (a:i--1) 



The first and second approximations to V^q/Vq for a few typical molecular models 
are given in the subjoined table : the model chosen is the n"' power centre of force, 
which is equivalent to a rigid elastic sphere when n = 00. 

Table III. 



Vo/V. 



Maxwellian 

gas, 



n 



5. 



1st approximation . 
2nd approximation . 



1-000 
1-000 



n 



9. 



1 • 000 
1 • 004 



n - 17. 



1-000 
1-008 




In calculating the approximations to V^^/V^ for n — 9 and n ~ 17 we have adopted 
the values 0'420 and 0*390 respectively for I2 {n)!!^ {n) ; it would be very troublesome 

■^ [^Added June ^, 1916. — It is interesting to notice that in all the cases in Table II. the error of the 
third approximation to Dq/D is equal (to the order of accuracy of one unit in the last figure given) to the 
difference between the second and third approximations. 

After each step in the approximation the ratio of successive increments seems to become progressively 
larger, and from Table I. it would be difficult to judge whether Dq/D converges to 1 * 12 or 1 * 13, if we did 
not know the correct value to be 1 - 1317. 
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to calculate the exact values of these functions, but fortunately the present calculation 
does not require them to be known with any great accuracy. When ti = 5 we have 
I2(5)/Ii(5) ~ 0*501, and when n — ^ the value is 0"333 (c/. Jeans' treatise, 2nd ed., 
§33), so that the assumed values when n =^ 9 and n~7 cannot be materially 
in error. 

The value of V^V^^ is in all the cases considered (ri = 5 to t?. i^ oo) very nearly 
equal to unity, the correction introduced by a second approximation being so small 
that further approximations are not likely to lead to any but a negligible increase in 
accuracy. Thus the exact value of V^^/V in the case of rigid elastic spherical 
molecules, for instance, is not likely to differ from 1*017 (slightly greater than 1*015) 
by more than one part in a thousand. 

The exact expression for Dn corresponding to molecules of the type just mentioned 
is consequently given by the following equation : — - 

(13*22) Dn = 1*017 ^ 



o 
'J 



2i/oo-^ {2h'wm)'' 



(13-23) = ^'^^^^ 



(13-24) . =1-200-, 

where in the last line we have made use of the formula for the coefficient of viscosity 
K for a simple gas, which has already been given by the author {' Phil. Trans.,' A, 
vol. 216, § 11 (D), p. 337), viz., 

(13*25) K '-=- 1*016 



In my first paper on the kinetic theory (' Phil. Trans.,' A, vol. 211, p. 477, 1912) 
the formulae (13*22)-(13*25) were given as above except for the omission of the 
factors 1*017 and 1*016 in the first and last, which resulted in 0*150 taking the place 
of 0*1520 in (13*23). 

The expression (13 '23) for Dj^ agrees almost exactly also with a result obtained by 
PiDDXJOK'^ for the same quantity, by an entirely difi'erent method. Mr. Pidduck's 
work is based on ITilbert's transformation of Boltzmann's integral equation for the 
velocity-distribution function. t His formula for Dn (^oc. cit., p. 101, 41) is the same 
as (13*23) except that the numerical constant, there given only to three places of 
decimals, is 0*151. 

* PiDDUCK,^ Proc. Lond, Math. Soc.,' (2), 15, p. 89, 1915. 

t HiLBERT, ' Math. Ann,/ 72, p. 662, 1912; "Boltzmann, ' Yorlesungen liber Gastheorie,' I. 
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We may also compare (13*24) with the result obtained from the simple mean-free- 
path theory (first introduced by Maxwell"^) with Jeans' correctionf for the 
persistence of velocities, viz. (in our notation) 

(13*26) Dn = 1*34 -. 

P 

This is therefore about 10 per cent, larger than the exact value (13*24). 

{g) The Variation of D^g with the Relative Concentration of the Component Gases, 

In the course of the development of the kinetic theory perhaps no branch has been 
the subject of more dispute than that dealing with diffusion. The point of greatest 
difference was the effect on the rate of diffusion of the relative concentration of the 
diffusing gases. Meyer's elementary mean-free-path theory J led to a formula for 
Di2 according to which the coefficient of diffusion should vary with the proportions of 
the mixture over the extreme range indicated by the equation 

(lS*27) \-^12/vi = __ ^2 

(Dl2)..= ^l' 



12 



Thus, when the molecular masses are very unequal, the range in the value of D 
should be very great. No such large variation is found to exist, however, according 
to the results of experiment. 

Meyer's theory took no account of the tendency of a molecule to persist in motion 
along its original direction after collision : as Jeans§ has shown in connection with 
viscosity, however, ^' persistence of velocities " is a very important fact, the neglect of 
which leads to grave error in the mean-free-path theory. Kuenen || has shown that 
when taken account of in the theory of diffusion, it largely removes the discrepancy 
between the small observed variations of D^g, and the variations theoretically 
calculated by the method referred to. 

An earlier modification of Meyer's theory by GROSsf may also be mentioned. This 
predicts variations of amou[nt similar to those observed, but generally of the wrong 
sign ; its merits are not such as to demand more than this brief historical reference. 

* Maxwell, * Scientific Papers,' !., p. 377, or 'Phil. Mag,,' 1860, January-July. 

t Jeans' ' Dynamical Theory of Gases,' p. 273. The whole of the chapter on diffusion (ch. xii. in the 
second edition) is of great interest, and a general reference may be made to it for comparison both of 
theory and experiment with the results of this memoir. 

I Meyer, ' Kinetic Theory of Gases ' (English edition), p. 255. 
§ Jeans' ' Dynamical Theory of Gases,' pp. 276, 292. 

II KuENEN, ' Supp. No. 8 to the Communications from the Leyden Physical Laboratory,' January, 1913. 
Of. also Jeans' treatise, ch. xiii. (2nd ed.), p. 328. 

H Gross, ' Wied. Ann.,' 40, p. 424, 1890; the disagreement of Gross's theory with experiment has 
been indicated by LoNiUS, 'Ann. d. Phys.,' 29, p. 664, 1909. 

2 B 2 
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In sharp distinction from Meyer's theory, Stefan'^' and Maxweli/I" put forward a 
theory of diffusion which led to the result that D^g is entirely independent of the 
ratio j^i : j^2? the formula involving only the sum V1 + P2 or p^. Stefan's theory was 
based on the hypothesis of rigid elastic spherical molecules, and agrees with the 
corresponding form of our first approximation to D12 (13'06) ; the same i*esult was 
subsequently obtained by Langevin| and the present author § on the most general 
molecular hypothesis, in the form ( 13*06). Maxwell's first formula related to the 
elastic sphere theory, and was greater than that of Stefan by one- third, Maxwell 
being considerably hi error : his second formula related to Maxwellian molecules, and 
was identical with the exact formula for such molecules given in § 13 (&). All these 
results share the property that they are independent of vi : ^2, and as the present theory 
shows, all of them are in error in this respect, except Maxwell's second formula. 
They require to be multiplied by the correcting factor V^/V of (13*03) : this factor is 
unity in the case of Maxwellian molecules, so that for a mixture of two sets of such 
molecules Dyj is independent of their numerical proportions. This is one of the few 
properties of a gas which depends in its very nature {i.e,, not merely in absolute 
magnitude) on a particular molecular model : another such property wilLbe noticed 
in § 14. The fact that the absence of variation is in disagreement with experimental 
results confirms the conclusion drawn from other sources of evidence that Maxwellian 
molecules are unsatisfactory representations of actual molecules. 

The only other general case of independence of D12 or vj : v^ is that in which the 
dynamical properties of the two sets of molecules are alike, as in the case of self- 
diffusion (§13 (/)). Hence we may speak of the coefficient of self-diffusion of a gas 
without specifying the ratio Pi : P2 in contemplation ; in general, on the contrary, Dy^ 
is defined only when the value of p^ : P2 is assigned. 

In the case when ms/m^ is negligible (§13 (e)) the mutual encounters of the light 
molecules are neglected, with the result that in this case also Dj2 depends only on p^y. 
This, however, is perhaps hardly to be regarded as an additional exception to the 
general rule. 

The error in the formulae of Stefan and the other authors mentioned arose from 
the neglect of the difference between the actual law of distribution of the peculiar 
velocities of the molecules and the assumed Maxwellian law. This is taken account 
of in the present paper, and leads to the correction factor V^JV in (13*03); we will 
now proceed to consider how this factor varies with the ratio p^ : P2 by studying the 
variation of the second approximation to it, i.e., (l— €o)~\ From this we may readily 

* Stefan, ' Wien. Siteb.,' 63 (2), p. 63, 1871 ; 65, p. 323, 1872. 

t Maxwell, ' Scientific Papers,' i., p. 392 ; ii., p. 57 and p. 345. Of, also Boltzmann, ^ Wien. Sitzb.,' 
66 (2), p. 324, 1872; 78, p. 733, 1878; 86, p. 63, 1882; 88, p. 835, 1883. Also ' Vorlesungen iiber 
Gastbeorie,' i., p. 96. 

i Langevin, loe. cit. ante, 

§ * Phil. Trans.,' A, toI. 211, p. 449. 
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infer the character of the variation of vyv itself, in view of our knowledge of the 
mode of convergence of the successive approximations as indicated in § 13 (e). 

(h) The Variation of (l— cq)"-^ ivith v^ : V2. 

So far as regards our detailed discussion of €o, we shall confine ourselves to the case 
of rigid elastic spherical molecules ; for Maxwellian molecules e^ is identically zero, 
while for n^^ power centres of force ( 5 < n < oo ) eo is intermediate between and the 
value appropriate to rigid elastic spheres {of. Tables II., III., § 13), the character of 
its variations being similar in the two cases. 

From § 9 (c) we have 



(13-28) eo = (^1-1)' 



^1^1^+ 2&i2^ll/2 +^2^2 



2' 



and in the case of rigid elastic spherical molecules {c/. § 9 (c) and § 9 (/)) we may 
write 



(13-29) 








^12 


(13-30) 


{h- 


-ly 




2 

Ml 


" 30 If > 17mi' 


where 










(13-31) 








P ~- M1M2 


so that 








' 



Q-A^p 



^ +78-H^p 



9 
P 



{k,-iy^ = 



J 2 

O9 f^2 



10.0^^1^.. 2 5 



d2 so -~~^^p- 17 JUL 






(13*32) 0<p<i 0<^<1. 

Since the suffix 1 refers to the heavier molecules, we have /tii > /xg and 

/1 q'QQ\ ^1 \^ ^2 

\ X t/ 00 J ~r~ ^-^ ~r- . 

(J/\ 0/2 

The condition that -—■ shall exceed -j- is readily proved to be 

Cii2 ^2 

(13-34) ^ j.^^^l-~-M2p-^2mp^ 

^ V /xi *^^^^ 4^+81p-180p'' 

and it is found that, for the admissible values of p {ie., 0<<p<i-), f{p) is positive 
and steadily diminishes as p increases; its least value is consequently /(|-), which is 
equal to 13/(9 + 4^). Since k never exceeds unity, /(i>)> 1, and hence M2i</(jp), so 
that 

(13*35) ^>^' 

(a/12 Oj2 
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The sign of equality corresponds to the limiting case k ~ 1, ^11=^ fj.^^ i.e., to the case 
of self-diffasion, in which instance, therefore, bjdi = 63/(^2 ~ ^is/^is? % being indepen- 
dent of vi : 1/25 as we have already seen (§ 13 {/)). 

We cannot unconditionally make the further statement that bifdi^by2/di2, since 
the necessary condition for this, viz., iu''i/jut2^f{p)^ is not always fulfill63d ; but it is 
easy to prove that 

(13-36) |i>|^ 

if 



12 



12 



{13-37) . /^>i^(l + 80M2-100MiM2). 

Ml 



When yoci = /X2 = J, the last condition becomes ^^1, so that if 7^; = 1 we have, as 
before, bjdi = bi^jdy^l if Mi = /^2 ™ i and ^ < 1, then b^jdi < byjdy^- The condition 
(13'37) may be roughly expressed by saying that the more equal the molecular 
masses, the more equal, also, must be the molecular radii in order that b^jd^ shall 
exceed &i2/<ii2 ; or, converse^, the more unequal the masses, the more unequal, also, 
may be the radii, consistently with the truth of (13' 37). 

Clearly, if (13*37) is satisfied, and 

di di2 d2 

€q steadily increases as the proportion of the heavier gas varies from to 1, and 
consequently, also, D^g steadily increases {cf. (13-07)). If, however, (13*37) is not 
satisfied, cq and D^g will first increase to a maximum and then diminish slightly, as 
the ratio v^ : vq increases from to 1. 

As regards the actual value of €q, the range (corresponding to all possible ratios of 
molecular mass or radius) is from the minimum value of {hi — lfcjf^ which is clearly 
zero (when /^cg = and vi ~ 0, i.e., when the heavier molecules are infinitesimal in 
number and the lighter molecules are infinitesimal in mass) to the maximum value of 
{k^ — iyb/e, which is ^^^3 (when /a2,^ and v^ ^ O). The maximum value of {hi — lfajd 
is 2 9^5 which lies between the above limits. The corresponding range of (l— £o)"\ the 
correction factor to the first approximation (13*06) to D12, which is introduced on 
making a second approximation, is consequently from 1 to 1*083. Hence, in con- 
junction with § 13 (e), we may conclude that for rigid elastic spherical molecules the 
total possible range in the complete correction factor V^jV to the first approximation 

to D12 (c/. (13*01)) is from 1 to ~ or 1*132 {cf. (1314)). 

It would not be difficult to construct a table showing the values of (l— eo)"*-^ for 
various typical pairs of gases, but owing to the fact that there are three variables 
concerned (^'.e., mjm2, ci/cra, and vjp2) it vv^ould need to be somewhat complicated, and 
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probably little more knowledge would be gained, from it than may be deduced from 
the above discussion and from Table lY. of § 13 (i), 

(i) Variation of D12 with j/^ : V2. Comparison ivith Experiment. 

In view of the conflicting theories regarding the relation between D^g and j/j/z/g, 
several series of experiments to test the question were made at Halle during the years 
1904-1909, under the direction of Prof E. DoRisr.'^' The gases used were helium — 
argon by Schmidt and Lonius, oxygen — hydrogen by Deutsch and Jaokmann, 
carbon dioxide — ^hydrogen and oxygen— nitrogen by Jagkmann. In Jackmann's 
summary of the experiments it is pointed out that the observed variations in D12, 
though far from negligible, are small compared with those predicted by Meyer's 
theory, while Gross's theory {cf. § 13 (g)) is equally unsuccessful in that it predicts 
variations of the right order but in the wrong direction. 

jEANst compares these observations both with Meyer's theory and with the 
corrected form of Meyer's theory after allowing for the persistence of velocities 
{cf. § 13 (g). The following table is taken from §446 of his treatise (2nd edition), 
and relates to the pair of gases, helium — argon : — 



Table IV. — -Variation of D12 for He- A. 



n/^2' 


D12 observed 
(Schmidt and Lonius). 


D12 (calculated). 


Meyer's theory. 


Meyer's theory corrected 

(Jeans, Kuenen). 


2 • 65 
1-00 
0-31 


0-961 
1-000 
1-036 


0-548 
1-000 
1-526 


0-910 
1-000 
1-110 



In this table the values of D12 have been multiplied by a factor in order to make the 
value corresponding to p^ = V2 equal to unity, for convenience of comparison. Jeans 
remarks, a propos of the above table, that the observed variations of D12 are 
insignificant compared with those predicted by Meyer's theory, being small even 
compared with those predicted by the corrected theory, and that Maxwell's simpler 
formula (according to which D12 is independent of 1/1/1/2) is after all the most accurate 
numerically. We shall see, however, that the present theory, with its correction 
V7V to the usual first approximation (13*06) to D^g, gives results which are in much 
closer accordance with the observed phenomena. 

^ ScHMmT, ^ Halle Dissertation' (1904), and 'Ann. d, Phys.,' 14, p. 801, 1904 j Deutsch, 'Halle,' 
1907; Jackmann, 'Halle,' 1906; Lonius, 'Halle,' 1909, and 'Ann. d. Phys.,' 29, p. 664, 1909, where 
the results of the whole series of experiments are summarised and discussed, 

t Jeans' ' Dynamical Theory of Gases ' (2nd ed.), § 446. 
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The following table contains the results of the Halle experiments on Djs, as given 

1 

in Jackmann's summary. Additional columns give the values of (a) -, calculated 

from (13*28) on the hypothesis that the molecules are rigid elastic spheres'^ : (6) 
VyV — -cf. ( 13*01)— estimated from the value of in accordance with § 9 (a) and 

1— eo 

§13 (e) : and (c) of D^g, obtained by multiplying the values of V'/V by such a factor 
as to make the mean of the resulting values of D^s equal to the mean of the 
corresponding set of observed values. This is all that is here necessary, since we are 
concerned now only with the variations of D^j, and not with its absolute values. 



Table V, — -Variation of Dio with 



12 



Vl I V2' 



Gases. 
Authority. 



Argon-helium. 



ScHMiDTf and LoNius, 
mi/w2 == 10. 



Oxygen-hydrogen . 

DEUTSCHt and Jackmann, 
milm^ =16. 



Carbon dioxide-hydrogen. 

Deutsch, 

mijrn^ ^ 22. 



Nitrogen-hydrogen. 



Jackmann, 
mil 1112 = 14. 



Oxygen-nitrogen . 



Jackmann, 
mim^ = 8/7. 



i 








, Di2. 


Ai - ""i 




1 


V' 






VI -\- P2 




1 ~€o 


v; 














Observed. 


Calculated. 
0-248 


f 0-273 




1-035 


1-050 


0-244 


0-315 




1-040 


1-056 


0-250 


0-250 


J 0-377t 
j 0-500t 




1-044 


1-063 


0-250t 


0-251 




1-050 


1-073 


0-254t 


0-254 


0-677 




1-058 


1-087 


0-256 


0-257 


[ 0-763 




1-062 


1-094 


0-263 


0-259 
j 0-276 


r 0-25t 
0-25 


j 


1-035 


1-050 < 


0-276t 
0-276 


^ 0-5t 
0-5 


■ 

1 


1-050 


1-073 1 


0-280t 
0-273 


j 0-282 


I 0-75t 




1-064 


1-097 


0-289t 


0-289 


1 0-25 




1-029 


1-040 


0-214 


0-212 


)> 0-5 




1-059 


1-089 


0-218 


0-222 


0-75 




1-070 


1-108 


0-228 


0-226 


1 0-235 




1-030 


1-042 


0-268 


0-263 


r 0-5 

i 

I 




1-050 


1-073 


0-266 


0-271 


1 0-467 




1-012 


1-013 


0-0733 


0-0731 


r 0-5 




1-013 


1-014 


0-0730 


0-0732 


-X 










1 



'^' The values of the molecular radii used in the calculation of €o were taken from the table on p. 476 of 
my first paper (' Phil. Trans.,' A, vol. 211, 1912); cf. the column there headed "elastic sphere theory." 
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We may conclude from the above table that the calculated values of Dig are in very 
satisfactory agreement with the observed data, the differences between the two being 
not greater than the experimental errors would render probable, if we may judge of 
these errors from the internal accordance of the observed results. It should also be 
remembered that the -molecular model chosen for calculation is not the best 
representation of an actual gas molecule, though it is sufficiently good for the 
purpose, especially in view of its simplicity for numerical work. 

(j) The Absolute Magnitude of Dig. 

We next consider the absolute magnitude of the coefficient of diffusion. A 
comparison of a theoretical expression for Dig with the corresponding observed value 
involves not only the accuracy of the theory and experiment, but also the suitability 
of the molecular model adopted as the basis of the theory. If we choose the rigid 
elastic sphere as model (and for many purposes this very simple model is fairly 
satisfactory) we may deduce from the observed value of Dig for a specified gas- 
mixture the corresponding value of o-i + crg (c/. (l3*0l) and (9*31)). By doing this for 
three pairs of gases A-B, B-C, C- A we can obviously determine from the resulting 
values of 0-^ + 0-5, o-^ + oTa, cTc + a-a the individual molecular radii cr^, o-^, cr^.^ By taking 
different sets of pairs we may obtain more than one determination of each molecular 

radius, and the mutual accordance of these affords some sort of check on the theory 

* 

— mainly, I think, relating to the suitability of the molecular model. But we cannot 
in this way test whether the theoretical formula is in error by a factor which is 
nearly or quite constant, since this would merely alter the deduced values of the radii 
in a common ratio. 

\_Revised June 2, 1916, — ^By another method, as follows, we can to some extent 
check the absolute magnitude of the theoretical results. Having determined values 
of a- in the above manner, we can use a-i (the radius for a particular gas) to calculate 
Dii by the formula which expresses Dig in the case of dissimilar molecules, putting 
mg = mi, (Tg = (Ti. In this way, practically by interpolation, we obtain a virtually 
experimental value of Dn (which cannot be measured directly). But the theoretical 

expression for Dn can also be written in the form J)\, = h^^, .. and p, being 

respectively the viscosity and density of the gas, while ^ is a numerical constant. If 
the correct theoretical expression for k^ is used here, the theoretical value of k should 
agree with the experimentally measured value Dnpi/zcp The former depends, of 
course, on the molecular hypothesis adopted, varying from 1120 for rigid elastic 
spheres to 1*504 for Maxwellian molecules. The following table gives several 
experimental values of h (taken from Jeans' treatise t), which all lie between the 

^ This method is due to Lord Kelvin, * Baltimore Lectures,' p. 295. 

t Of, Jeans' ' Dynamical Theory of Gases,' 2nd ed., §§ 447, 448. The value of h deduced from the 
corrected Meyer's theory is 1*34 as against the Talue h =^ 1-20 gmn by the present theory (13-24). 
The formula for k there used is th^t giyen in mj second memoir. 
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values here theoretically deduced. The values of n, calculated from the variation of 
viscosity with temperature, on the hypothesis that the molecules- are the n^^ power 
centres of force, are given alongside. The formula of interpolation used to obtain D^ 
from Di2 is based on the theory of rigid elastic spherical molecules, and this, together 
with the rather early date of the diffusion constants used (they were determined by 
Loschmidt), may prevent us from expecting a very close agreement with theory. I 
have not calculated the theoretical values of k for values of n between 5 and oo^ 
since they depend on certain integrals which must be evaluated by quadratures ; but 
it would be interesting to have these quantities determined.] 



n. 



DiiPi 



Ki 



Theory . . . , 
Hydrogen . . . 
Carbon-monoxide . 
Oxygen .... 
Carbon-dioxide 
Theory . . . . 



« • 



00 

12 
9-3 

7-9 
5-2 
5*0 



1-200 



1 
1 



36 
34 
1-43 
1-50 
1-504 



Perhaps the most satisfactory comparison between theory and experiment is 
provided by calculating a-i both from the coefficient of diffusion, as above, and also 
in some other way as, e.g., from the coefficient of viscosity /ci. Assuming that the 
experimental data and the molecular representation underlying the theory are 
satisfactory, the accordance of the two values of cr affords a complete check on the 
theory of D^g and k jointly. Such a comparison was made in my first memoir (' Phil 
Trans.,' A, vol 211, p, 480), using the approximate formulae for D^2 and k there 
given ; the confirmation of the theory thus provided was perhaps as good as the 
accuracy of the data might allow one to expect. The table in that paper requires 
revision, because of the corrections to my original formulae which have been deduced 
in the present and my second memoir. Until more modern and accurate data for 
the coefficient of diffusion are available, however (and of such there is need), this 
revision is hardly worth while in view of the smallness of the corrections mentioned. 

In connection with this matter, finally, it may be pointed out that in tables giving 
values of Djg intended for accurate theoretical purposes it is. necessary to specify the 
ratio by volumes of the diffusing gases, since the theoretical expression for Byj, in 
general, involves this ratio (fj : V2). This practice has not usually been adopted 
hitherto. 

(k) The Variation of D12 with Pressure and Temperature, 

Since the theoretical expression for D^g, on whatever molecular hypothesis, contains 
at least one quantity (molecular diameter or force constant) which is not directly 
measurable, two or mora values of D^g, corresponding to different temperatures or 
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pressures, are necessary in order to obtain evidence from diffusion as to the best 
molecular model. It is, therefore, of some importance to consider the variation of 
Di2 with pressure (or density) and temperature. 

The first approximation to D^g (13'06) varies as T2(w-i)+^ in the case of molecules 
which are n^^ power centres of force, the case n = ^, T'^% corresponding to rigid 
elastic spheres (c/. ' Phil. Trans.,' A, vol. 211, p. 479). It may readily be seen that 
the correction factor VyV in (l3'0l) is independent of A or T for these particular 
molecular models {cf, ' Phil. Trans.,' A, vol. 216, § 9 (B), p. 321). 

Another molecular model, considered in my two previous memoirs, is the rigid 
elastic sphere surrounded by a field of attractive force (the Stjtherlani) molecule). 
The first approximation to D^g in this varies as T'^'/(l4-S7T), where S^ is known. as 
Sutherland's constant of diffusion."^ It may be proved without difficulty, as in the 
case of the coefficient of viscosity for a simple gas {cf, my second memoir, §§ 9-11) that 
the correction Vyv is usually intermediate between unity and the value corresponding 
to rigid elastic spheres without attraction. The correction factor in the case of 
Sutherland's molecules also depends upon the temperature, but the variation is so 
very slight as to be negligible. 

As regards variation with density, the first approximation to D^s varies inversely 
as the total density of the gas mixture, and is independent of the relative proportions 
of the two gases. The correction factor VyV in (13*01 ) has been seen to vary with 
the relative proportions {vijv^^ but if vi : v^ is fixed, the factor may readily be shown 
not to vary with the total density v,-¥v,. 

§ 14. The Coefficient op Thermal Diffusion. 

{a) The General Formula for Drp. 

The coefficient of thermal diffusion Drp was defined in § 10 (a) by the equation 
/ \ / -r-x 1 o J_ * "i~x (J lo£r JL 

(1401) t6A= — JJrn — t;;^ — ^^^ '-' -L'T Ti 5 

1l ox ax 

and was found to be given by {cf, (lO'lO), (5'36)) 

(14*02) Urp = ---9l>0Poi = -9X1X21301 -i^rr^ j T* 

In the case of Maxwellian molecules it was found that ^^o vanishes {cf, (6*05)), so 
that for a gas composed of such molecules D,j is zero, and the phenomenon of thermal 
diffusion is non-existent. This, together with the absence of variation of I)i2 with 
nh2 {cf- § 13 {g)), is one of the few properties of a gas which depends essentially on 

* Cf. Sutherland, ' PMl. Mag.' (5), 38, p. 1, 1894, and my first memoir, p. 479. 

A \j Zi 
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the particular characteristics of the molecules. We shall see later that the 
phenomenon is also difierential, i.e., it occurs only when the masses or diameters of 
the molecules are not the same for the two kinds. Thus when they are alike in mass 
and dynamical properties /3^o (and consequently also Drj) is exactly zero {cf, (8M3)) ; 
this is true whatever the nature of the molecules. 

(b) The First Approximation to D^. 

The first approximation to D^ is obtained by substituting the value of /3^o given, to 
this degree of accuracy, by (9-20), viz. (of, (917), (9-09)-(913)), 

(14-03) J),= Hrn^-^rnjm ejk^^ 

27n/omim2K 12 (0) AiVq 

We shall not go beyond a first approximation in the general case, but (as in § 13 (e) 
with regard to D^a) we shall briefly examine the order of magnitude of the corrections 
introduced by further approximations, by a consideration of the special case when 
m^/mi and 0-2/^1 are very small (§ 7). 

(c) The Value of Dfj; when mi/mg and 0-1/0-2 are very large. 
In this instance the value of /3'o is given by (7 '22), and 



1 



(14-04) D, = -^^|^g. 



When the molecules are rigid elastic spheres of radii o-j, o-g, this assumes the form 



(14-05) J)^ = 



3X2 J)\ 



16vo{cr,~\-cT2y{h7rm2y^'I)' 



or (of. (7*26)) the same expression except that BJD replaces DVD^ In each case the 
suffix 1 to D or D' denotes the second minor of the corresponding determinant. 

We have already given the value of D as far as the fourth row and column (§ 13 (e)). 
From this we find that the first three approximations to BJD are, with their 
successive differences equal to 

2280 ^ 
4608 

30 
556416 



1061681 



0-524. 
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These indicate the magnitude of D^/D sufficiently for our purpose, though for an 
accuracy of (say) J per cent, it might perhaps be necessary to carry the calculation 
two or three steps further. If we take the limit of these approximations as 0'58, 
we have (c/ (l3*13)) . 



(14-06) 



Dt - 0-58 



3X, 



16i^o {^i + ^2Y {hirm^^^ 



0-51XoD 



2 



12- 



In this case, therefore, D^j is comparable with D^g, though it may be remarked that 
if X2 is too nearly equal to unity, our assumption (§7) that the effect of collisions 
among the lighter molecules is negligible may require revision. It may, indeed, 
be readily seen from the general expression for D^ that this always vanishes if either 
Xj or Xg becomes zero. 

We may briefly examine the case also of n*^ power centres of force by means 
of (13*19). Considering only a first approximation, the following results are obtained 
for various values of n : — • 



n — 


5. 


(m : 


1) 


n — 


7 1 


'm ■ 


-f) 


n ~ 
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(m : 
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13 


(m 
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17 


(m 


+) 


n — 


33 


(m 


= -k) 


n 


00 


(m 


-0) 



First Approximations to T)\IW, 



J)\/D^ = 0, exactly. 



0*156, approximately. 

0*227 

0*294 

0*326 

0*372 

0*417 



The general formula for the first approximation to D^^/D^ is readily seen to be 



5 
2 



■m 



6—m 



12 



U' 



n- 



5, 
3 



/ 
/ 

V 



m 



n- 



The last table shows that a very slight excess of n over its value for Maxwellian 
molecules (^ = 5) suffices to raise Dr^ to a considerable proportion of its magnitude 
corresponding to ^ = co (rigid elastic spheres). The phenomenon of thermal diffusion 
clearly disappears only under conditions which must be fulfilled with great nicety. 

As an instance of the correction introduced by a second approximation in the case 
of a finite value of n^ it may be noticed that the result of a second approximation 
for n == 9 is 0*246. 
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(d) The Approximate Value of D^ when a-i — (Tg and m^lm2 is very large. 

The case considered in the last sub-section is realised physically in a metal or gas, 
among the molecules of which are found free electrons. Another case of physical 
importance is that in which the molecules are similar except in mass, mjm2 being 
large ; the diameters of ordinary molecules do not differ greatly, while their masses 
may be widely different. Thus the diameters of the oxygen and hydrogen molecules 
are approximately 178 . 10"^ and 1*33 . 10"^ respectively, the ratio of their masses 
being 16. If in our first approximation to Drp given by (14*03) we insert the values 
of K^s (O), ^1, &c., appropriate to rigid elastic spherical molecules, and write o-^ ~= o-^, 
/xi ~ 1, /;t2 = 0, the result is 

^' = 4{9 + X.)l6.„{., + .,r \h^ ' approximately. 
The corresponding value of eo(915) is 3Xi/(39 + X2), so that 

-n 39 + X2 3 fm.-\-mX^' . ' -, 

^^^ = W^) 16.0 (C.. + ..)- \h^^J ' Woximately, 



an( 



(14'07) Dt = ~~-^— ^Di2? approximately. 

Thus when Xi = X2 = -g-, Dt/Di2 becomes 0*095, which is nearly equal to the 
maximum value for this particular gas. When Xi or X^ vanishes, so also does D^j., 
another instance of the general theorem on which remark has already been made. 

(e) General Remarks and Numerical Values. 

It is now evident, from the special cases treated in the preceding sub- sections, that 
the importance of the phenomenon of thermal diffusion, as compared with ordinary 
diffusion, is greater or less according as the difference between the molecular masses 
and diameters is large or small. If the mass and diameter of the one set of molecules 
are very large compared with those for the other set, D^ rises to about I-D13 ; if there 
is a large inequality in mass, but no difference in size, D^ is about i^yDjs at most^^' ; 
if mass and diameter are both the same for the different molecules, Drj. ~- 0, which 
is also its limiting value when X^ or X2 tends to zero, and its invariable value when 
the molecules are Maxwellian. 

The following are a few values of Drj./Dj,2 corresponding to typical pairs of actual 



•^- 



See Note E (p. 197) for the case of nearly equal masses and very unequal sizes. 
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gases, and are obtained by comparing the approximations to D^ and D12 contained 
in (13'07) and (14"03); calling their ratio Ic^, we have 



k 



T 



gl (l— g()) 



V, 



0^1 



2 (^1 - 1 ) [1 5 (X1/X2 + X2M1) (mj — M2) — 4miM2^i2^ (^1—^2) + {\^i2ku -^^21^2)] 

fM]M2 (^] —^2) + 6 (Xa^/— ^W) + (^2Ml-A-i— XiM2-A.2) 

AA + (6-^3ViM2) (A, + A,)— i^(mMi+m2'A2)+-^/-WmiM2 ' 



where 



£!l^ 



4 _Zl 
3 



/ 20-1 






2a-2 



1/2/X2 Xo-i + 0-2/ 



-2—3 



^iMl Vi + O-gy 



The values of o- used in calculating the results in the following table are taken from 
the table on p. 476 of my first memoir. 

Table VI. — Comparison of D^, Dp, Dig. 





Argon-helium. 


Oxygen-hydrogen, 


Oxygen-nitrogen . 




^2 


™i 16, 

^2 


^1 _ 16 


Eatio of the 


^1:= 1-69, 

org 


^1 = 1-34, 

0^2 


^1-0-97, 

0-2 


two gaseous com- 
ponents. 


7 20-4Xi + 4-8X2 


7^ 28-0Ai-f-5-2X2 


j^ __ 0-33Xi4-0-4U2 




32-6 + 58-6-V + 4-4-2 


A/rjl - •• 

45-54- 78-3''i +5-8'^^ 


22-2-1-9-3-1 +8-7-^2 




V2 1^1 


z/2 ^'l 


1^2 ^1 


k - ^- 


k - ^P 


Dp 

Dt* 


'^^"■Di2- 




Dp 

Dt' 


-L'12 


^ Di2- 


Dp 

Dt* 


A, - 1 ^1-1 

s/2 


0-133 


0-513 


3-9 


0-122 


0-592 


4-8 


0-008 


0-026 


3-2 


^1 23 1 

1/2 


0-132 


0-408 


3-1 


0-128 


0-441 


3-4 


0-010 


0-033 


3-3 


Ai-I, "i-f 


0-079 


0-166 


2-1 


0-079 


0-230 


2-9 


0-007 


0-024 


3-4 


^2 





















We may remark, finally, on the direction of thermal diffusion : since D^ has been 
found to be a positive quantity (I4'0l) shows that the heavier gas tends to difiuse 
towards regions of lower temperature, while the lighter gas tends to diffuse towards 
those of higher temperature.^ 

^ Of,, however, Note E, p. 197. 
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§15. The Coefficient of Pressure Diffusion. 

In § 10 we saw that, among other causes, a pressure gradient is an agent in 
producing diffusion, according to the law {cf, ( 10*08) 



(15-01) «; = D,^^ = D,^l2gPo 

^Pq ax ^ ax 

where (c/ (10*09) ), 

(15-02) D^ = "^ B,, = ^^^A^i-^^) D^^ ^ L D 



'p-^m 



if 

(15-03) . ^ X,X,(m,-~m,) ^ 

The phenomenon depends essentially on the difference of mass of the molecules, and 
we may note also that, like thermal diffusion, it tends to zero with A^ or X^. If mi/ma 
is very large, k^ clearly has the value X2? except when X2 is very nearly equal to 1. 

It is interesting to compare the degree of diffusion produced by equal gradients of 
log ^0 and log T. This is equivalent to comparing k^ and krj^. When milni2 is very 
large, kjkrj^ — l/0*51 or approximately^ 2, if (T^lcr2 is also large {cf. (14'06)), or 
(39+X2)/l5Xi if cTi — (T2 {cf, (14*07)), save when Xi is small. The factor k^, unlike /c^, 
is not dependent on the relative size of the molecules. 

In Table VI. are given the values of kp corresponding to the three typical pairs 
of gases considered in § 14 (c). It appears that D^ is of thrice or four times the 
magnitude of D^ for such gases. 



§ 16. The Steady State Without Diffusion. 

We will now briefly consider the steady state, without diffusion, of a gas subject 
to the influence of {a) external forces or {h) boundary conditions which maintain a 
constant non-uniform distribution of temperature. If we write zero in place of u^^ in 
(10*11), and divide by D12, we have 

(16-01) ^^-2h(T,+ ^^)+h^-^ == 

OX \ Vq ox j dx 

as the equation of state. It is more convenient to express the middle term as in the 

following equations (substituting for - -^^ from 2*04- — — being zero, since the state 

vq ox at 
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of the gas is steady — and expressing the result in terms of Xj and X2 instead of Xo 

(16-02) ^ = 2hxM^i-^2)-h^^^^. 

(16-03) 1^^ = 2h\X2{X,^X,) + kr, ^-^^^ 

Since 

(16'04) 



9 i^iS 3 /XA 1 3Xt Xi 3X2 _ XJ+X2 3Xi _ 1 3x, 






dx\v2/ dx\X2/ X2 3a^ Xg^ 3cc Xg^ dx X^ '^x^ 

we may re-write (16*02) or (16*03) in the form 

(16-06) 3 log (.-,/'.) ^ 24 (X.-X,)- -^ ^J^. 

d^ XjX2 d^ 

If the temperature and external forces are uniform, from (16*06) and the two 
similar y and z equations we may deduce the integral 

(16*07) "^ = p) 6'M<3:i-X.)x+(Y.-Y,V+(Zi~-Z.)4^ 

V2 \V2/0 

which is a well-known result* directly deducible from the statistical theory of 
distribution of density in a gas."^ If the external forces are gravitational, the z axis 
being yertically upwards, 

= 1=0, Zi = —niig^ Zi2 = —'^^2^, 

and, if g be regarded as constant, we have 

V2 \vJo 

The heavier gas, naturally, is relatively denser in the lower strata, the amount of the 
effect being greater the greater the inequality of mass and the smaller the temperature. 
In the case of the atmosphere, since the molecular masses of oxygen and nitrogen are 
nearly equal, the magnitude of this imperfection of diffusion is but small. It is found 
that the change in the value of 1/1/1/2 would amount only to about ^ per cent, per 
kilometre, t 

"^ Of. Jeans' * Dynamical Theory of Gases ' (2nd. ed.), p. 91 (234, 235). 
t Cf. JiUNS' treatise, § 369. 

VOL. GCXVII. — -A. 2 D 
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The variation of relative density in the atmosphere due to the temperature 
gradient (which is about 7° 0, per kilometre) is still smaller, and in the opposite 
direction. Thus, taking vi : v^ for oxygen and nitrogen as 1 : 3 roughly, and T as 
280° C. approximately, the equation giving the variation of v^ : v^ near the earth's 
surface (independently of the effect of the difference in mass) is 



8 lv\_ 1, 13T 

3% 



dx\vol x/ T9x' 



3 



By Table VI., /% — O'OOB, X^ = f, while if we adopt one kilometre as tlie unit of 
length, writing T = 280 and ™~- as —7, the change in - per unit length is found 

to be 

— • ^"008 . -4- • 7 = 0-00035, 

9 280 

or ^ per cent, per kilometre. This variation is of one-sixth the magnitude of the 
gravity variation, which it tends to neutralise, since the heavier gas tends to rise 
owing to the temperature gradient. These changes are masked, however, by the 
effects of convection currents ; it has been found that the proportions of nitrogen and 
oxygen in our atmosphere are the same for all altitudes up to 10 kilometres, 
and, indeed, above this height, right into the isothermal layer or stratosphere. 

As another instance of the phenomenon of a concentration gradient due "to 
temperature, we may consider the case of a mixture of equal parts (by volume) of 
oxygen and hydrogen, or argon and helium, placed between two plates which are 
maintained at freezing and boiling points respectively. Neglecting the small 
variation in v^jv^, as far as the integration with respect to x is concerned, in the 
equation 

we may conclude that, independently of the distance between the plates,"^ the change 
in the value of i^i/j/g is 

~-4.0-130.1og;f^| ^ 0161, 

where we have taken h^ --=^ 0'130 {cf. Table VI.), which is approximately the case 
either for Ar-He or 0-H. Thus near the hot plate there will be approximately 
48 per cent, of the heavier gas to 52 per cent, of the lighter, and mce versd near 
the cool plate. This remarkable result is independent of the pressure of the gas, 
assuming this to be uniform between the two plates, as will be the case in the 

* To the order of accuracy according to which hj. may be regarded as constant. 
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absence of external forces ; the density, on the other hand, will vary inversely as the 
temperature. 

These numbers are based on the hypothesis that the molecules are rigid elastic 
spheres ; for n*^ power centres of force the effect would be somewhat less, vanishing 
altogether if n ^ 5, It seems very desirable that the theory should be experimen- 
tally tested ; the effect predicted is of easily measurable amount, and could be further 
magnified by taking a greater temperature range. The magnitude of the concen- 
tration gradient is so large as to seem improbable, and it is possible that some 
circumstance has been overlooked which would modify the theory, but I have been 
unable to detect any such flaw. It is difficult to say how long a time it would take 
to reach the steady state, or what influence the constant flow of heat through the 
gas, from the hot plate to the cold, would have upon the phenomenon. (See Note D, 
p. 196.) 

§ 17. The Coefficient of Viscosity. 
The general expression for the coefficient of viscosity is (c/. (ir06)) 

(17*01) Kyj '■=^ ^|-5- -— 2 (l^irr + ^2y-r)' 



e shall not trouble here to go beyond a first approximation to this expression, 
even this being somewhat complicated. Referring to (9*27), we obtain the following 
result as our first approximation : — 

SRT Civi -4- 2c^2nv'i ■+■ <^2^2 



(17*02) 



^12 



ttK'is (0) FnCii/i' 4- 2h\^G\^v^v^ -f h\^C^v^^ ' 



The values of Ci, Cg, Ci2, c\2 are defined by (9*28)-~(9'30). The equation (17 '02) is 
identical, except as regards notation, with that given in my first memoir on the 
kinetic theory {loc, cit,, p. 451). When V2 ■=■ 0, the formula reduces to 

^,j,.Qg\ 5ET 1 5RT 



xK\,(0)V" ^KnMo)' 

which is the first approximation to the coefficient of viscosity of a simple gas 
composed entirely of molecules of the first kind. Some idea of the accuracy of 
(17'02) as a first approximation may be gained from the fact that for a simple gas 
the error (i.e., of (17'03)) amounts to only 1'6 per cent, for rigid elastic spherical 
molecules, and less for m}'^ power centres of force (c/. my second memoir, loc. cit., § 11), 
the first approximation being too small. How the error varies with v^ : V2 cannot be 
stated without carrying the calculation further, but it is probable that it is always 
of the same order of magnitude, 1 or 2 per cent. A second approximation to k^^ 
would replace the last factor of (17*02) by the quotient of one homogeneous quartic 

2 I) 2 
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polynomial in v^, v^ by another, reducing, when 1/3 = 0, to the second approximation to 
ici, the coefficient of viscosity of a simple gas. The present importance of such a 
second approximation to K12 does not seem to render the complicated calculation worth 
while. In my first memoir {loc. cit,, pp. 469-472) a comparison of (17*02) with 
experimental data is given, which suggests that the formula satisfactorily represents 
the behaviour of actual gases in respect of the variation of /c^^ with pjvs.^ 
When the molecules are rigid elastic spheres, (17 '02) has the special form 

-pm 1 ir (1+1^12) Vl^ + i \ 1 + -Z P~7~^ \ ^1^2 + (1 +i^2l) ^2 

li 1 m{m^ \^^ t 2^i/x2 (0-1 + 0-2) J _ 



where we have quoted from § 9 (/). The special forms appropriate to molecules of 
other particular types may likewise easily be written down. 

§ 18. The Coefficient op Theemal Conduction. 

In § 12 the following expression for the coefficient of thermal conduction 3- was 
found (c/ (12-21)):— 

(18-01) ^ - i-B,^ ji(^,/3, + .2^-.r)-^ 2r(.,a, + .2fit_,)|. 

oil a 1 J 

Owing to the complexity of even a first approximation, we shall not go further 
than this in the present paper ; in making the approximation, we use the results 
of § 9, which lead to the equation 

(18-02) 3- :^ ^^ (^1^1 + ^2^2) (^1 + ^2) C.RT 

27rmim2K^i2(0) ip^i'^Y -^'^(^12^1^2'^ ^2^2^ 

M2 2,0^ , Ml 2\ 2 \kx—l) VxV^^l 

L\yfc.n V / di<4-2rfi2^i^2 + rf2VJ 

If we put vx = in this expression it becomes 



(18-03) 9- - .-^^J^^^-^ =: I^^C, 



JilTfCii Jx 12 (0/ 



by (17'03) ; this agrees with the first approximation to 3- for a simple gas which was 
given in my two previous memoirs {loc. cit,, p. 462, p. 337, respectively). The 

* In particular, the theory there p.Yen indicates that the viscosity of a gas mixture may rise to a 
maximum (for a certain value of vifv^ which exceeds the viscosity of either component separately. This 
had already been proved by Maxwell for the particular type of gas dealt with in his second paper. 
KtJENBN has recently shown that the ordinary elastic-sphere theory leads to a similar result, when 
allowance is made for the persistence of velocities after collision ('Amsterdam Acad. Proc.,^ 1%^ p. 1162, 
1914). 
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expression in (18'02) for a composite gas does not, however, agree with the corre- 
sponding first approximation in my original memoir (^ Phil. Trans.,' A, vol. 211, p. 452, 
cf. 40 and 4l). This is because it was there assumed (§2) that thermal conduction 
was independent of dijffusion, whereas we have seen that a temperature gradient 
is inseparable from either diffusion or a concentration gradient ; this invalidates the 
equations (24) of that paper. The second term within the square bracket of (18*02) 
especially arises from the terms in the velocity- distribution function which represent 
the diffusion effect , although the first term is also affected. Evidently even the first 
approximation to % involves the quotient of a homogeneous quartic in v^, v^ by the 
product of two homogeneous quadratics, a complexity which makes it hardly worth 
while to compare the result with the scanty experimental data at present available. 
A brief discussion of my earlier formula, which may be expected to differ only slightly 
from (18'02) in numerical magnitude, is given in § 18 of my first memoir. 

5 19. The Specifio Energy of Diffusion. 

In § 12 (c) we saw that in a gas which is at rest and at a uniform temperature, 
so that no conduction of heat is taking place, there may yet be a continual rise in 
temperature if diffusion is going oil The gas being at rest as a whole, no gain 
of thermal energy accrues through the medium of internal friction, and the thermal 
flux of diffusion may proceed also in the absence of external forces, ^.e., solely as 
a result of variations of concentration. However the diffusion is produced, in the 
latter or any of the other ways described in § 1 0, the equation of energy will contain 
a term equivalent to that on the right of the following equation :— 

(19-01) - ; poa^-s^(iiO. 

This is identical with (l2*22), and assumes that only the thermal flux of diffusion 
is operating to increase the heat energy of the gas. 

By (12*23) we have the following expression for H, which we have termed the 
specific energy of diffusion : — 

(19'02) © = ™~^— Sr(i/ia,+F2a_,) 

« J <X 1 



J/0 RT Drj 



X1X2 J Dj2 

by (9*20) and (10*05), (lO'lO). Since (X, the specific heat at constant volume for the 

composite gas, is given by 

R 



(19-03) ^''^^Jm 
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we have, writing hr^ for Drji/Di2, as in § 14, 

(19*04) ^^^^SH^^vl^^ 



A1A2 

Thus we may write (l9'0l) in the form 
(19-05) p,G^^ = _ 2 A || _|. ^^C„T< 

We may notice that there is an interesting analogy between this equation and the 
ordinary equation of continuity 

.(19-06) 1^ = - 2 ^ (potto), 

which refers the rate of concentration of matter in a volume element to the differential 
rate of transfer of matter across its extent. In the case of diffusion, the opposing 
inter -diffusing streams carry equal numbers of molecules in opposite directions, but 
while the energy of molecular agitation is the same for either group of molecules, the 
kinetic energy of diffusion is different, owing to the different masses with the same 
velocity of diffusion ± {u^q, v^q, w^q). It is, I imagine, to this cause that the thermal 
flux of diffusion is due, though the suggestion is only tentative and does not affect 
the accuracy of (19*05) one way or the other. It is clear, however, that if the 
velocity of diffusion is such as to cause a concentration of the more massive molecules 
in any region, at the expense of the lighter molecules (since the total number v^ is not 
affected by diffusion), the temperature of the region will rise, ^^ being positive (§ 14) ; 
this readily follows from, a comparison of (19*05) and (19'06). From our discussion 
of the magnitude of kj^ in § 14, it appears that the specific energy of diffusion is 
greater the greater the difference in mass and size between the molecules, as we 
should expect. 

§ 20. Appendix on the Inequality oe Temperature between the Component 

Gases. 

Added June 2, 1916. 

Up to the end of § 3 the equations of this paper take account of the possibility of 
a difference of temperature (as defined in § 1 (c)) between the component gases. 
After that point it is virtually assumed that T^^ is zero. We will now briefly 
consider what modification must be made in order to cover the general case. 

It is clear from (3*151) that the component gases have different temperatures only 
when the ratio of the two gases, by volume, is changing with time. Consequently in 
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all steady states the temperatures are the same, although the gas may be non-uniform 
in temperature, velocity, or composition ; in this case §§ 4-19 give a complete account 
of the first-order phenomena of the gas. 

When the relative proportions of the component gas are varying, so that the 
temperatures are unequal, the determination of the coefficients a, ^, y in §§ 4-9 is 
unaffected, but we must examine how far the physical equations of the succeeding 
sections remain valid. It is immediately evident that all the phenomena grouped 
imder the general ternl ''diffusion" {cf, § 10 and §§ 13-16) are uninfluenced by the 
presence or absence of the temperature difference, because the former depend solely 
on the terms in f/{f)o which are of odd degree in U, V, W {cf. § 3 (a)), while the 

terms relating to T'o ^^ "t~^ ^^^ ^^ ^^^^ degree. 



Further, it is readily manifest that the equations of viscosity are likewise unmodified. 
The mean hydrostatic pressures of the component gases pi, p^ are indeed altered, but 
their sum p^ remains constant with the value defined in terms of v^ and To by (I'lS). 
Since the series J{C^) of § 3 (a) is symmetrical in U, V, W, the partial pressures 
{Pxx)i} {Pyy)i^ {Pzzji ^^e all affected equally with pi, and similarly for the corresponding 
pressures for the second component gas. The differences p^x-^Pm Pyy~^Pii^ Pxz-^Po are 
therefore independent of T^o, and it is clear from (3*03), (3'04) that this is true also 
of jpy^, p^^., p^y. Hence the equations of viscous stress given in § 11 are applicable 
both to steady and unsteady states of the gas. 

Finally, it may be seen on inspection that the equation of energy, deduced in § 12, 
is true independently of the existence of a temperature difference between the 
component gases. The equations (12*06), (12*07), on which the equation of energy is 
based, remain true in all cases, if T is taken to be Tq, Also the presence of the even 
power series J (C^) in the expression for /(U, V, W) does not affect the value of UC^, 
so that the expressions of § 12 (c) for the coefficient of thermal conduction and the 
specific energy of diffusion are universally true. 

It remains only to form an estimate of the magnitude of T^o by determining 
approximately the coefficients §, in J(C^). For this purpose it is simplest to make 
use of equation (2*23), viz., 






2s 



Then it may readily be proved, after the manner of §§ 6, 7 of my second paper {loc. 
cit), and as a consequence of the expressions (3*03), (3*04) for /(U^ V^ W), that 

(20-02) - ^ ^f'^^J—rx ^^i"^ == 32Do|^^ 2 [{dn^^{r,s) + d,Ur.s)} Sr+d,n2{r.s) S_^,l 
1.3.5... (2,9 + 1) at r = o 

(20*03) ^—.i^!^^^ = 32Do^ i [d,,,,{r,s)§^+ {d,,,,{r,s) + d,,,,{r,s)} S^^l 

^ 1.3.5... (2^+1) at r = o 
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where 
(20*04) 



dim {r, s) 



2 



r-l,s~in 



r, fi 



g-(.^+y') ^2y2 2 ^jj* (y) B\,,, (v, s) dx dy, 



k~0 



M J» 



(20-05) d,,,, {r, s) = N,_,, ,_^,v,v, e-<^'+^'> x^ 2 <f,,/ (y) B\,,, {r, s) dx dy, 



r, S 



(20-06) 4n.(n s) - N,_i.,,^i.i.3 e-'^'^y^xY 2 {-\Y<i>r}{ti)W,,,,{T,s)dxdy. 



r, s 

2 

^=0 



There are similar equations in which the suffixes 1, 2 are interchanged, which need 
not be written down. 

From a comparison of (20°01) with (20*02), (20°03) we obtain the following 
equations for 4? ^-^r •~™- 



(20*07) 



(20*08) 



1 ^ 
32Do- 2 [{dun{r, s) + di,r2{r, s)} <5, + (^2n2(n s)S_,'] = 1, 

Vo 

•i 00 

32Do- ^[d,22i{r, s)S,+ {d2222{r, s)-{-d2i2i{r, s)} S_,] = -1. 

Vo 



In these equations s ranges from 1 to oo (that the zero value is excluded may be 
readily seen from § 2 (c), (e)). Similarly r effectively ranges over the same values, 
since it is clear from (20*04)-(20*06) that d{0, s) = in all cases.'^ Hence the 
equations (20*07), (20*08) do not enable us to determine 4 and S___q, which are given 
by (3*101), (3*102) in terms of the remaining §'s. Indeed, to obtain T^o '^^ terms of 

-— ^, by means of (3*151), we need only to calculate Zr^ or E rJ_^ {cf. (3*122)). 

Symmetrical determinantal expressions for these can be deduced as in § 5, if desired, 
but we shall be content here to determine S^q to a first approximation only. 

The two central terms of the two central equations of type (20*07), (20*08) are 



(20*09) 



(20*10) 



1 



32Do-[{c?nn(l, l) + d,2u{l. I)}^i + rf.n2(l, l) S^,] - 1 

^0 






S2J),~-d,22i{l, 1)^1+ {32222(1, l) + 4i2i(l. l)}^-~i] = -1, 



and it may readily be seen that 



(20*11) 



(20*12) 



dini{hl)=7^.^i'K\,iO) 



144 



^2222 V 1 J 1 / 



TT 
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l^2^K^22(0), 



(20*13) (^1212(1? 1) — — C?i22l(lj 1) = <^212l(l? 1) '■= •~'^2112(l> 1) = 5^ /^lM2^1^2K^12 (0 

00 



* When r or s is zero, h == 0, and f/>^'^(y) ^-^ 0. (7/1 (4*17). 
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Hence, if we assign to Dq (which has so far been arbitrary) the value given by 

(20-14) Do=- \, 

we may write the equations (20*09), (20*10) as follows : — 

(2015) (l + ^Vi-^.i = 1, 

(2016) -§,+ (l + ^^)§.^ = -l. 
We have here introduced the notation 

/9nn ^\ Z- = -—jyii^/ Z. = ^ 22 (Q) 

^"^^^^^ "' -K'i,(0)' ' -K',,(0)- 

From (3122) we deduce the following equation for (5'o : — 



00 



(2018) ^\= -2\,\^'Zr{8,-^_,). 

1 

The first approximation to this is found, from (2015), (2016), to be given by 

(2019) ^'o = -2XiX ^12^1'-^ ^sA' 



2 h^h^* 

^12«^l +^21 "^2 + 



4miM2 
Hence, by (3151) and (2014) we have 



(20-20) T'o - iD-To^^o 1^ 



dt 



3To [j^ h%' Y'd\\ 



47n/oMlM2K^12 (0) I 4^iM2 (^12^1° + >^21^2^) J 9^ 

IT) ^1+^2 [-| , ^l'^2' I"' 9X^0 

^ '' 2R I 4MxM2(XiA' + X,,yfc/)J dt 



I 



3X^0 



where Di2^ is the first approximation to the coefiicient of diffusion (c/. (13*08)), while 
\ is a new quantity, defined by the last equation, which we may term the 
anisothermal difiusion constant. 

We may easily gain some idea of the order of magnitude of I^ if we consider two 
gases of similar dynamical properties, Le., such that m^ = mg, and k^ =. h^ ~ I. The 
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last factor in I;^ is nearly equal to unity, so that we shall neglect it in our approximate 
calculations. In this case, therefore. 



T — 1 
1a = — ^ 



D 







For oxygen the value of Dj/ is approximately 0*19 {of, Jeans' treatise, p. 335), 
while 'R/m = 260.10^ {loc. ciL, p. 131). Hence I;v = 0*24 . 10~^ approximately. 
Thus if the density-ratio of the two gases is changing by 1 per cent, per second, the 
difference of temperature, if they are equally mixed {i.e,^ Xi = Xg = J), is approxi- 
mately one thousand-millionth of a degree centigrade."^ The hotter gas is the 
one which is diminishing in relative concentration. In fact, we may perhaps regard 
the slightly excess speed of the molecules of this gas as the cause of its relative 
expansion and consequent rarefaction. It is clear from these figures, however, that 
the phenomenon is likely to be of very minute dimensions in ordinarj^ cases. 

^ ThusT'o = XiA2(Ti~T2) = ^(Ti-Ta) - -0-24.10-7^^ = - 0-24. 10"^ so that Ti -Tg = - IQ-^ 



U 



approximately. 



Note A, (To p. 118, line 3.) 

[Since this paper was communicated to the Royal Society, this part of the theory has been 
experimentally tested and confirmed, at least qualitatively. An account of the preliminary qualitative 
experiments, made by Dr. F. W. Dootson, will be found in the ' Philosophical Magazine,' xxxiii, 

p. MS.— Note added Fehruary 22, 1917.] 

Note B. (To p. 118, middle.) 

[Recently I have succeeded in proving the identity of tt and sin rrx with the determinants referred 
to, using elementary algebraical methods. The construction of a formal justification of the analytical 
methods used in the present paper is also being proceeded with, and I hope to be able later to deal 
with the questions of convergence raised here, and to bring the present theory into satisfactory relation 
with that based on Boltzmann's equation. — Note added Fehruary 22, 1917.] 

Note a (To p. 129, line 2.) 

[In a later paper it will be shown that considerations of convergence require a re-grouping of the terms 
in these expansions, so that they become series of polynomials in C^. The expansions are really used in 
this form in the present paper, the re-grouping of the terms being effected by means of the difference- 
transformations described in § 5 (d). — Note added Fehruary 22, 1917.] 

Note D. (Top. 189, line 12.) 

[Since the above was written the indicated result has been qualitatively confirmed in the case of three 
pairs of gases. The time taken to attain the steady state was only a few hours in the experiments 
referred to ; they are described in a Note by Dr. F. W. Dootson and the author in the ' Philosophical 
Magazine/ xxxiii, p* 248.— iVo^^^ added Fehruary 22, 1917.] 
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Note E, (To pp. 158, 184, 185.) 

[In a note recently communicated to the ' PhilosopMcal Magazine,' I have considered thermal diffusion 
in the case of two molecules nearly or quite equal in mass, and very unequal in size. If the difference of 
mass is sufficiently small, the larger molecules will tend towards the cooler regions. — Ajpril SO, 1917.] 

NoteF. (Top. 115.) 

[In his 'Inaugural Dissertation,' Upsala, 1917 (received just before the final revision of the proofs of 
this paper), D. Enskog gives a mathematical theory of simple and composite monatomic gases, based upon 
Boltzmann's integral equation for the velocity-distribution function. . The method of solution is, however, 
different from that of his 1912 paper. The numerical and other results, including those relating to 
thermal diffusion, are in agreement with those of this paper, though not always identical in form. While 
perhaps less developed from the physical standpoint. Dr. Enskog's work is mathematically much the more 
complete. His elegant and accurate proofs will materially lighten the task of proving my own work to 
be in conformity with Boltzmann's equation. — A;pril SO, 1917.] 



